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ABSTRACT  (ONGERUBRICEERD) 

A  rniall  lidar  system  has  been  used  to  measure  tibe  vertical  structure  of  the  atmosphenc  extinctioa  in  a  dune 
area  bordering  The  Hague  (The  Netherlands),  at  about  2.6  km  from  the  North  Sea.  The  atmospheric  optical 
properties  at  diis  location  are  detennined  by  a  mixture  of  industrial,  urban,  rural  and  marine  aerosols,  which 
composition  depends  on  die  air  mass  history.  The  measurements  were  made  unattended,  around  the  clock, 
five  days  a  week.  About  250  extinctimi  profiles  were  recorded  every  day.  This  rqxnt  reviews  die  data  base 
obtained  and  presents  some  selected  results.  Tlw  hdar  ^stem  is  described  briefly.  Factors  influencing  the 
accuracy  of  the  inversion  of  lidar  signals  are  discussed. 
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SAMENVATTING  (ONGERUBRICEERD) 

Een  klein  lidaisysteem  is  ingezet  voor  het  meten  vaa  de  vertikale  atmosfeariscbe  extinctie  piofielen  in  de 
dumen  aan  de  rand  van  Den  Haag  (Nederland)  op  een  afstand  van  ca.  2,6  km  van  de  Noordzee.  De  optische 
eigenscbappen  van  de  atmosfeer  aldaar  wordm  bqwald  door  een  combinatie  van  industriele,  stedelijke, 
natourlijk  continentale  en  marine  aerosolra,  afhankelijk  van  de  histone  van  de  hicbtmassa.  Het  gehele  proces 
van  meten,  data  opslag  en  data  verweridng  vond  crmtinu  en  automatisch  plaats  gedureade  vijf  dagen  per 
week.  Per  dag  werden  ongeveer  2S0  extinctie  profielen  gemeten.  Dit  rapport  geeft  een  overzicht  van  de 
verkregoi  gegevens.  Een  selectie  van  de  resoltaten  wordt  gqnesenteerd.  Het  gehniikte  lidarsysteem  zal  in 
het  kart  wordoi  beschreven.  Factoien  die  de  nanwkeurigheid  b^vloeden  van  de  inversie  van  lidarsignalen, 
wotden  besproken. 
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INTRODUCTION 


In  the  period  January  1985  to  August  1986,  vertical  extinction  profiles  of  the  x!lmoq)liere  were 
measured  with  an  automatic  lidar  system  (optical  radar)  at  the  TNO  Fbysks  oul  Electronics 
Laboratory  (52'.07.43"  N;  04M9'.40'*  E).  The  location  is  indicated  in  Hgure  1.1.  The  system  was 
mounted  in  a  cabin  on  the  roof  of  die  institute,  about  15  m  above  ground  level  and  poinling  in  a 
northern  direction  over  the  dunes.  In  princ^le,  an  extinction  profile  was  measured  every  6 
minutes,  5  days  a  week.  Each  {»ofile  measurement  was  preceded  by  one  or  mme  horizontal 
measurements,  to  determine  the  horizontal  extinction-  and  backscatter  coefficients. 
Meteorological  parameters  were  recorded  in  combination  widi  the  lidar  data.  Cloud  cover  and 
prec^itation  were  recorded  during  working  hours  (when  possible).  Qoad  base  height  was 
inverted  fiom  the  lidar  data.  Results  were  presented  in  Kunz,  1990.  In  die  first  few  mondis  of  the 
project  also  a  sodar  (acoustic  radar)  was  operated  but  the  results  were  unsadsfectory  to  continue 
these  measurements. 


Rgurel.l;  Sdiematicniqi  of  the  re^on.  The  measuremem  site  is  indicated  by  a'*'. 

Large  indastrial  areas  are  dotted. 
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The  lidar  opented  aotomatically  under  cooixol  of  an  HP>8S  personal  ccm^uter.  The 
meteorological  station  was  operated  widt  a  PET-4000  Ccmimodare  Business  Ccmqmter.  Bodi 
systems  were  programmed  in  BASIC.  The  lidar  signals  wero  processed  to  crfxain  the  horizontal 
backscatter-  and  extinctmn  coefBcknts  (at  a  level  of  IS  m  altitude)  and  to  obtain  tibe  vertical 
extinction  pn^les  using  die  forward  integratioo  mediod  (Kunz.  1983),  see  also  Chapter  3.  The 
inverted  results  were  stored  »«««»»■<<  of  die  raw  data.  The  main  reason  for  diis  was  die  limited  data 
storage  ciqiacity  at  diat  tune,  hi  a  later  stage,  die  data  was  converted  to  IBM-PC  format  for 
analysis.  Afterward,  the  lidar  and  meteorological  data  were  merged  into  a  data  base  wbii:fi  finaUy 
coaqnised  about  23,000  records.  This  data  base  provides  a  unique  source  for  correlation  of  lidar- 
measuied  extinction,  backscattw  and  vertical  profiles  widi  die  actual  meteorological  ccmdidons  to 
determine  statistical  relations  between  dieae  parameters.  Rud^rmore,  die  padi-integrated 
extinction  can  be  determined  as  a  function  of  elevation  angle.  Slant-padi  atmospheric  extinction 
plays  an  important  role  in  models  for  predicting  die  performance  of  electro-optic  sensors. 

This  report  describes  the  hdar  system,  its  calibration  and  the  inversion  techniques  ^lied.  Facuns 
are  discussed  which  influence  the  accuracy  of  die  calculated  extinction-  and  backscatter 
coefficients  at  low  altitudes.  An  overview  of  the  acquired  data  base  and  some  selected  results  are 
presented  and  discussed. 
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2  UDAR  SYSTEM 

The  TNO  Mini-Lidar,  diown  in  Hgine  2.1.  is  a  modified  laser  range  finder,  developed  by  the 
(Datch)  Oldelft  Optical  Industries  Qvesently  named:  Delft  Electro  0|:«ics)  as  a  military  laser 
range  finder  at  die  end  of  the  sixties.  The  powerful  laser  (100  mJ)  and  large  telescope  (74  mm 
diameter)  made  the  system  attractive  for  modification  to  a  lidar  system.  The  sy^em  qqieared  to 
be  very  reliaUe  and  has  been  in  operation  by  TNO-FEL  since  1982.  It  was  used  both  for  the  semi- 
condnuous  measurements  during  die  2>year  VISA  project  described  in  diis  report,  as  well  as 
during  several  national  and  international  field  experiments  on  land  and  over  sea.  The  heart  of  the 
Mioi-Lidar  is  a  NdrGlass  laser  with  a  rotating  prism  as  (^switch.  The  specificadons  of  the  lidar 
are  sumrrunized  in  Table  2.1. 


Rgiiur2.1;  1740  Mmi  Lidar 
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Table  2.1:  Pn^wities  of  die  TNO  Mini-Lidar 


Transmitter. 

wavelength 

1060  nm 

energy 

SOmJ 

pulse  duration 

30  ns 

beam  divergence 

Imrad 

repetition  rate 

6ppm 

Receiver: 

diameter  telescope 

74  mm 

focal  length 

6S5tnm 

RCA  C30916  APD  (diam.) 

13  mm 

sensitivity  @  1060  nm 

12A/W 

coIot  filter  RG  830 

3nim 

sensitivity 

19.3  kV/W  (500) 

noise  level 

0.16mV(— 8.310-’W) 

Laser  Energy 

The  energy  of  the  laser  was  measured  with  three  different  energy  meters.  The  calibration 
differences  between  die  meters  were  less  dian  2  %.  The  variation  of  the  laser  energy,  over  200 
shots,  was  less  than  3.3  %.  It  is  assumed  that  the  influence  of  the  energy  variation  over  the 
measuiement  periods  is  negligible.  During  the  project,  the  laser  energy  was  measured  daily  and 
used  to  calibrate  die  backscatter  measurements. 


Sensitivity  of  die  Receiver 

The  accuracy  of  die  receiver  has  been  tested  with  a  calibrated  light  emitting  diode  (LED,  1060  nm 
and  with  a  power  level  of  194  nW;  Kunz,  1988  and  1989).  The  emitted  light  was  collimated  to  fill 
die  field  of  view  of  die  receiver. 


RmggftfthpLidar 

hi  general,  the  range  of  a  lidar  depends  on  die  properties  of  die  system  (e.g.  laser  energy,  receiver 
diameter  and  noise  of  the  d^ector)  and  on  the  actual  weather  conditum.  For  our  Mini-Lidar,  die 
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theoretical  maximum  range  is  about  1.3  km,  in  a  homogeneous  atmos{diete  with  an  extinction  of 
0.4  km'^  In  a  molecular  (Rayleigh)  atmoqdtere  die  maximum  range  is  about  4S0  m.  The  method 
for  calculating  die  maximum  range  of  lidar  systems  in  a  homogeneous  atmo^ete,  is  mitlined  in 
Appendix  A.  In  cases  of  an  inhomogeneous  atmosphere,  the  range  dqiends  on  die  structure  and 
cannot  be  loedicted  a  priori. 


Beam  Direction 

The  slant  path  direction  of  the  beam  in  the  Mini'Lidar  can  be  controlled  via  the  internal  mirror, 
from  an  angle  of  about  -10*  to  about  -<-33*  with  respect  to  the  horizon.  To  reach  a  maximum 
altitude  with  the  system,  an  external  mirror  under  45*  was  qiplied  during  the  first  few  months  of 
the  experiments.  However,  during  periods  of  prec^tadon  this  resulted  in  erroneous  signals 
caused  by  scattering  at  the  droplets  on  die  minor.  Therefore  it  was  decided  to  apply  only  the 
internal  minor  of  the  system,  giving  up  the  altitude  advantage  for  all-weather  operations. 


Pre-processing 

A  logahthmic  anqilifier  (Analog  Modules  LA-90'P)  has  been  used  to  suppress  the  dynamic  range 
of  the  measured  signals.  This  unit  has  a  logaridimic  transfer  over  the  first  80  dB  of  the  input 
signal  strength  (0  dB  =  5  V)  and  a  constant  transfer  for  smaller  input  signals  ( U,<  5. 10~*  V).  The 
results  of  the  calibration  and  some  odier  properties  of  this  unit  are  described  in  Appendix  B. 


Dm  Recording  and  Processing 

The  transient  recorder  was  deseed  and  constructed  in-house  at  the  beginning  of  die  project 
because  at  that  time  no  instruments  were  crunmerdally  available  with  the  desired  accuracy, 
resolution  and  sanqile  rate.  The  resolution  of  die  transient  recorder  was  9  bits  and  die  san^le  rate 
was  16  MHz.  The  trigger  accuracy  was  better  than  8  ns  to  inqnove  die  timing  accuracy.  Errors  in 
the  calculated  backscatter-  and  extinction  coefficients,  caused  by  small  time  shifts  of  the  lidar 
signal,  are  described  in  Appoidix  C.  An  HP-8S  personal  computer  controlted  die  lidar  system  and 
processed  all  data. 


7NOr«port 


Page 

12 


3  INVERSION  METHOD 

3.1  buioductioii 

This  paragiq>h  describes  in  short  die  model  of  die  lidar  signal  and  the  methods  used  to  invert  die 
signals.  A  conqnehensive  discussion  <m  this  subject  can  be  found  in  the  literature  (e.g.,  Klen 
1983,  Kunz  1983).  Some  qiecial  topics  are  discussed  mote  extensively  because  they  have  a  large 
impact  on  the  accuracy  of  the  results  and  have  not  been  described  elsewhere. 

The  general  model  for  lidar  signals  is  given  by  the  following  equation; 
where: 


Pd(R)=E, 


c  3(R)  A  _ 
2'  4ic  R^ 


R>0 


A 


El 

Pd(R) 

R 

m) 

T 

^oja 

OKR) 


area  of  die  leceivo'  in  km^ 
velocity  of  light  in  km/s 
lasCT  energy  in  J 
lidar  letum  on  the  detecttv  in  W 
range  in  km 

backscatter  coefficient  in  km'^ 
transmission  of  the  cities 
spatial  extinction  coefficiem  in  km'^ 


(3.1) 


The  duratitm  of  the  laser  pulse  does  not  appear  in  diis  fonnula  but  must  be  short  in  comparison 
with  the  smallest  range  to  be  observed.  The  receiver  and  the  transmitter  field  of  views  are 
assumed  to  have  conqilete  overliqi. 


The  ouqmt  of  the  receiver  in  volts  is: 

Ui{R)  =  DPa(R) 
where: 


(3.2) 
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Ui(R)  =  output  of  the  receiver  in  V 

D  =  sensitiviQr  of  the  jdioto  detector  in  V/W 


Throughout  the  project,  the  lidar  signals  were  pre-processed  with  an  analog  logaridunic  anqilifier 
to  compress  the  large  dynamic  range  of  the  signals.  Thus,  the  signals  could  be  recorded  with  a 
transient  reccrda  with  a  smaller  dynamic  range  (at  the  cost  of  the  reliability  of  very  low 
extinction  values).  The  prc^terties  of  the  logarithmic  anqtlifier  are  described  in  Appendix  B.  The 
output  signal  of  the  logaridunic  anqtlifia  can  be  described  by: 


U„(R)  =  Li+L2-ln[Ui(R)/Uj 


(3.3) 


where: 


Li 

Uo(R) 

Ur 


intercept  of  the  transfer  in  V 
slope  of  die  transfer  in  V 
ou^t  voltage  of  the  log-anqilifier  in  V 
reference  voltage  of  1  V 


Division  by  U,  is  only  a  formal  action  to  render  the  argumoit  of  the  log  fiincticui  dimensionless. 


Lidar  signals  can  be  divided  in  two  classes,  i.e.  firom  a  homogeneous  atmosphoe  and  from  an 
inhomogeneous  atmosphere.  It  is  assumed  that  the  lowest  layer  in  the  atooosidiere  is  well  mixed 
(over  the  range  of  the  lidar)  which  results  in  a  homogeneous  distribution  of  the  aerosols.  When 
the  atmosphere  is  sensed  in  the  vertical  direction,  however,  the  structure  will  goieiaUy  be 
inhomogeneous. 


The  purpose  of  lidar  is  to  determine  the  atmo^heric  badcscatter-  and  extincdon  coefBcimts.  In 
principle,  these  two  quantities  cannot  sinqily  be  determined  with  cme  single  ^stenL  However,  in 
a  homogeneous  situation  it  is  possible  to  determine  diese  two  quantities  because  the  strength  of 
the  lidar  signal  is  a  measure  for  die  backscatter  and  the  time  behaviour  gives  die  extiru^on  (see 
section  3.2).  fri  situations  of  an  inhomoyeoeons  atmos[dieie  the  ariqilitude  and  die  time  behaviour 
of  the  lidar  signal  ate  determined  by  die  spatial  properties  of  both  die  extinction  and  the 
backscatter  coefficients.  In  diat  case,  additional  information  (a  relation  between  extinction  and 
badtscatter  as  well  as  a  boundary  condititm)  is  required  to  invert  the  signal. 
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The  techniques  used  for  the  inversion  of  Udar  signals,  eifber  from  homogeneous  or  from 
inhomogenecMis  atiho^dieies,  are  discussed  in  sections  3^  and  3.3.  Parameters  that  effect  die 
accuracy  of  die  results  are  also  discussed.  Because  die  effects  are  easiest  to  understand  in  a 
homogeneous  atmosphere  most  of  diis  discussiMi  ajipears  in  section  32.  The  results  from  the 
inhomogeneous  atmo8|diere  are  also  sensitive  to  these  errors,  but  their  effects  can  mily  be 
quantified  when  infonnadon  on  the  shape  of  the  profile  is  available. 

3.2  Inversion  of  horizontally  measured  signals 

Lidar  signals  from  a  homogeneous  atmosjdiere,  widi  a(R>-a  and  in  equation  (3.1),  are 

first  conected  fm  die  range  term  in  the  denominator.  Taking  die  logaridim  of  die  range 
conqiensated  waveform  yields: 


ln[S(R)]  =  In  (K,)  +  In  0)  -  2.a.R  (3.4) 

The  sl(^  and  the  intercept  of  equation  (3.4)  are  determined  by  linear  regression  which  provides  a 
direct  measure  for  the  extinction  coefficient  (oc±fra)  and  the  backscatter  coefficient  iTi(prHAp) 
This  procedure  is  generally  known  as  the  'slt^  method'. 

The  standard  deviadrm  in  the  extinction  and  the  backscatter  coefficimts  must  be  critically 
considered  for  the  following  reasons.  In  a  very  dear  atmosphero,  the  signal  is  veiy  weak  and  thus 
the  range  is  limited.  Moreover,  die  slope  of  die  range  compensated  signal  is  very  small.  This  may 
lead  to  unexpected  large  uncertainties  in  die  calculated  extinction-  (od:Ap)  and  backscatter 
coefficients  [ln(p±AP)  ]n(P)  +  In(l-i-Ap/p)].  On  the  other  hand,  if  the  extixiction  is  large,  the 
signal  is  strong  but  the  range  of  the  Udar  signal  (and  dius  die  number  of  sanqiles)  is  also  liininvt 
due  to  attenuation.  (A  mediod  to  calculate  die  maximum  range  of  a  lidar  ^stem  is  presmted  in 
Appendix  A.)  Due  to  die  limited  number  of  samples,  diis  can  also  lead  to  large  uncertainties  in  the 
calculated  extinction-  and  badcscatter  t  efficients.  In  conclusion,  die  errors  are  smallest  for 
intermediate  extinctions,  hi  Rgure  3.1  die  standard  deviatitms  in  die  extinction,  Act,  in  die 
badcscatter,  (AP),  as  inverted  from  die  data  base,  are  presented  as  a  function  of  the  extinction  and 
the  backscatter,  reqiectivety.  As  diown,  die  deviation  of  die  largest  extinction  values  (10  km**)  is 
less  than  1  knr*  (Le.  10  %),  while  for  die  lowest  extinction  values  (0.05  knr‘)  diey  may  amoiitit  to 
10-100  %!  Best  results  are  obtained  widi  this  system  for  extinction  confidents  around  03-1  km** 
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which  is  about  a  factor  of  two  larger  dtan  du  extinctko  at  udiich  the  masimuin  range  is  obtained 
The  standard  deviation  for  die  badcscatter  shows  a  mitiiiniim  around  0.02  km'>  and  increases  with 
the  badcscatter.  A  more  detailed  analysis  of  die  errors  in  die  extinction-  and  badcscatter 
coefficients  detennined  the  slope  mdbod  is  presented  in  Knnz  and  DeLeeuw,  1992. 


LOO(EXTINCTION  in  km  ’)  LOG(BACKSCATTER  in  km  ') 


Figure  3.1:  Deviatkin  in  the  exdnctioa  coefficient  and  in  the  backscalter  cOeffickot  as  a  fiincticm  of  the 
values  the  extinction  and  die  badcscatter  ooeffidents,  during  die  first  months  of  1985.  Tbe 
maximum  lange  of  the  Mini-Lidar  is  about  1.3  km.  This  is  reached  when  the  extinction  is 
about  0.4  km*^  and  die  badcscatter  is  about  0.028  km'V  See  also  Appendix  A 

Aldiou^  the  slope  mediod  seems  to  be  very  sinqile,  many  effects  can  influence  die  final 
results.  The  resulting  errors  are  not  only  dqiendent  on  the  accuracy  of  the  measurement,  but 
also  on  die  strengdi  of  die  signal  itself.  Some  sources  of  enrors  are: 


Range  compensation  is  die  bask  operation  in  lidar  signal  processing.  However  udieo  die 
trigger  mmiient  is  not  accurately  known,  unintended  conapensatkn  with  (R±dR)  ^  is 
perfoimed  instead  of  widi  R^.  This  leads  to  errors  in  bodi  the  extincxion  and  the  badcscatter 
coefficients.  Fbr  mtaii^Ie,  an  error  of  one  sanqile  interval  of  30  ns  introduces  an  error  in  die 
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extinction  is  0.1  km'^.  The  efExa  of  die  trigger  momeDt  accuracy  are  discussed  in  moie 
detail  in  appendix  C. 

-  die  acoiracv  of  tlw»  tr«n«fcr  nf  rtw;  lf>y-«mnHfigr 

Tbe  transfer  fimctkia  of  die  logaridunic  amplifier  has  two  constants,  L|  and  L^.  The 

HnewrtaintW  in  the  enwrtan**  mfhiewee  hoA  the  herlnratter  and  die  extinction.  The  amw  in 
not  («ily  dqiendent  on  die  emv  in  die  transfer  but  also  on  die  stiengdi  of  die  signal  In 
particular,  low  extinction  valnes  ate  sensitive  to  nnceriaintiea  in  the  transfer  of  the  log* 
amplifier.  As  explained  in  Appendix  D,  die  uncertainty  in  die  backscatter  dne  to  the  transfer 
function  is  less  dian  3  %,  but  die  error  in  die  extinction  varies  from  less  dian  2  %,  for  1 
km-i,  to  mote  dian  300  %,  for  ors  0.01  knri. 

die  influence  of  the  bandwidth 

Bandwiddi  reduction  has  the  advantage  to  improve  die  signal-to*noise  ratio.  However,  die 
sh^  of  die  signal  is  distwted  when  die  bandwidth  reduction  is  bt^ond  the  frequency 
contents  of  die  lidar  signal  For  exanqile,  signals  widi  a  short  rise  time,  such  as  frmn  a  cloud 
return,  ate  reduced  in  amplitude.  Also  signals  from  a  homogeneous  atmosphere  are 
influenced.  Although  deccmvolutuHi  can  dieoretically  restore  the  original  signal  eiqierience 
leams  diat  large  errors  remain  in  die  final  results.  The  amount  of  error  introduced  dqiends 
both  on  die  strmgth  of  die  signal  and  cm  die  bandwidth  of  die  filter,  fr  is  not  sinqile  to  give  a 
closed  fmm  eiqnessum  for  the  error  because  die  actual  error  dqiends  also  on  the  cross-over 
function  of  die  tidar  and  on  the  initial  condition  of  die  filter.  Simulations  of  homogeneous 
signals,  with  extinctions  from  lO-^knr^  to  10km*i  and  widi  filter  time  constants  between  10 
and  1000  ns,  show  diat  a  range  of  about  2.  lO^FmCT  is  tbe  time  constant  in  seconds  of  the 
first  order  filter)  is  reguned  to  eliminate  die  influence  of  the  reduced  bandwiddi  (Kunz, 
1979). 

-  tbe  accuracy  of  die  reference  level 

Signals  are  tenqwtarily  stored  in  a  transient  recorder  which  must  have  die  cqiabili^  to  store 
small  negative  signals  as  well  (noise).  Hus  means  diat  die  rignals  are  recorded  widi  a  small 
(numerical)  offset  wfaidi  must  be  subtracted  before  die  data  are  processed.  When  diis 
subtraction  is  not  done  or  not  done  accurately,  large  errors  in  die  cakulated  results  ate 
possible.  Also  here,  a  closed  form  eiqHession  for  die  error  is  not  possible  because  it  depends 
Ml  die  system  in  questkm,  the  transfer  die  dectronics  and  die  actual  valnes  of  die 


extiiictMMi-mddielMdac«ttercoefficiaitK.TIiM|«>blewilM«Mlini»ii<Bii^u^«tHrfliy 

Shimnni  (1985). 


die  reiohitioD  of  A/n  converter 

An  inbereitt  propecty  of  lidar  signals  is  Ibe  laige  dynamic  noge  vdiidk  varies  from  3  to  6 
decades  or  more.  To  record  diese  signals  widi  sufficient  accuracy,  an  A/D  converter  of  10 
bits  or  more  is  required.  When  the  VISA  experiments  were  made  (1S185-1S186)  diese 
instruments  were  very  expensive,  if  available.  Part  of  the  problem  has  bem  solved 
iqiplying  a  logaridunic  anqilifier  widi  a  dynamic  range  in  the  order  of  43  decades.  It  has 
been  shown  (Kunz,  1990)  that  a  combinatinn  of  a  Ingaridimic  anylifigr  anrf  an  right  hit  A/n 
ccmverter  is  ctMnparahle  widi  a  fourteen  bits  A/D  converter.  Most  errors  occur  in  «mall 
signals  at  large  ranges. 

3.3  Inversion  of  vertically  measured  signals 

Inversion  of  the  extinction  profile  from  die  measured  data  is  possible  widi  e.g.  equation  (33) 
vhich  has  been  described  by  a  number  of  audiors  e.g.  HitschfeM  and  Bordan  (1954),  Klett  (1983) 
and  Kunz  (1983). 


a(R)=- 


S(R) 

IT 


where  S(R)  is  die  range-corrected  signal  as  shown  in  (3.4) 


(3.5) 


As  can  be  seen  from  (33),  a(R)  can  only  be  inverted  if  both  S(R)  a«t  a  boundary  condidon 
Q((Rq)  are  available.  Fuidiermore,  die  equation  can  be  solved  by  integratiog  from  die  fra- 
end  (R<R^  or  from  near-by  (R>R^.  The  forward  integrating  mediod  has  die  advantage  that  die 
necessary  boundary  condition,  o(R^,  can  be  obtained  widi  the  lidar  itself  from  a  measurement  in 
die  horizontal  direction  (where  the  atmoqphere  is  assumed  to  be  homogeneous,  hence  die  slope 
mediod  can  be  qiplied  to  die  measurement).  However,  dus  mediod  hu  die  disadvantage  that  a 
singulaiity  may  occur  when  die  boundary  conditioo  is  too  large.  The  backward-integration'  | 

mediod  has  die  advantage  that  diexe  win  be  no  singularity  in  die  solution  becanse  a  real  irqwt  | 


bomidny  of  the  extmctkm  is  ahn^  positive.  However,  s  disadvantage  of  die  backwaed* 
integratkn  method  is  that  the  detenniiiation  of  die  boaadaiy  vahw  from  the  lidar  aigaal  haelf  is 
not  accmate.  For  the  mvcrskn  of  lidar  as  described  in  this  iqiort,  the  frvwanl  inteytinn 

method  has  been  applied. 

Other  methods  to  measnte  the  vertical  extinction  and  badrsratter  prrtBes  am  described  in  Kma, 
1988  and  in  Knnz.  1990. 

The  effects  described  in  die  ptevkms  sectioti  also  introdnce  errors  in  die  extinction  profile 
inverted  from  lidar  measmements  in  an  mhomcigeiieoiis  atmoqihete.  However,  because  diere  me 
so  many  difierent  profile  shqies,  it  is  not  possible  to  even  indicate  die  effect  of  possible  errors  or 
to  describe  a  general  trend.  A  detailed  analysis  is  necessary  for  eadi  type  of  profile. 


4 


I^CRlPnON  OF  THE  DATA  BASE 


This  dispier  gives  n  overview  of  the  active  periods  of  Ae  bdar  system  and  die  number  of 
profiles  recorded  during  VISA.  The  lidar  data  and  the  indq[»eodendy  recorded  meteo  data  were 
first  validated  and  subsequently  merged.  The  parameters  recorded  in  the  final  data  base  are 
described  in  ^ipendix  E. 

Hgure  4.1  shows  a  cakadar  of  the  years  1985  and  1986  in  wfaidi  die  active  di^s  of  die  lidar  are 
indicated.  The  lidar  was  not  operated  during  holidi^s,  weekends,  field  eiqierimeots, 
reconstnictioo  activities,  data  cooverskm  and  maintenance  of  die  equipmeat  The  di^s  for  which 
bodi  lidar  and  meteo  data  are  available  have  been  encircled  in  die  calendar.  As  can  be  seen  fnnn 
Figure  4.1,  winter  and  qiting  periods  are  emphasised  as  regards  die  number  of  measurements. 
Very  few  measurements  are  available  in  summer  and  no  measurements  at  all  in  Sqitember. 
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Figure  4.1: 
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The  foUo'wing  tid>le  gives  i  muiMarioi  overview  of  tihe  number  of  messuianeots: 


Udar 

meteo 

year 

fiks 

rBOotds 

files 

reooidt 

1983 

118 

25,234 

103 

7395 

1986 

71 

12,231 

64 

3.110 

total 

189 

37,483 

167 

12303 

Note:  die  meteo  data  was  avenged  every  15  minutes;  die  lidar  measured  a  profile  about  every  7 
minutes  which  causes  die  discrqiancy  between  die  number  of  lidar  records  and  die  number  of 
meteo  records. 

After  validation,  merging  and  integration,  the  following  number  of  files  and  records  are  available 
for  analysis: 


Combiaed  files 

year 

files 

records 

1985 

91 

13,869 

1986 

46 

6333 

total 

137 

22,804 

Analysis  of  the  data  base  is  described  in  die  next  section. 
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5  INVERTED  DATA 

Some  typical  examples  of  horizontal  extmctkm-  and  badfscatter  coefSdeats,  in  relation  with  die 
actual  meteorological  condition,  as  well  as  some  results  of  vextkal  atmoqdieric  extinction 
profiles,  are  presented  in  dus  dacgtax.  The  altitude  of  the  dood  base  as  measured  during  VISA  has 
been  described  in  Kunz,  1990. 

S.l  Analysis  of  die  horizontal  measuremrots 

For  the  inversion  of  lidar  signals  to  vertical  atmo^heric  extinction  profiles,  die  horizontal 
extincdon-  and  backscatler  coefficients  ate  necessary  input  parameters.  They  are  also  measured 
widi  die  lidar.  A  logical  sequence  is  dierefore  to  concentrate  on  diese  parameters  first. 

distribution  of  extinction  and  backscatter 

Statistical  overviews  of  die  horizontal  extinction-  and  badcscatter  coefficirots  are  presented  in 
Hgure  S.l.  All  available  records  (22,804)  were  processed  and  distributed  over  200  channels. 


to*  10’  10*  10’  10*  10*  10*  10’  10*  10’ 

EXTINCTION  In  knr’  BACKSCATTER  In  km  ’ 

RgureS.l:  Fiequcocy  of  occurrence  of  the  22.804  lidar  measured  horizontal  exincticn 

(a)  and  backacaiter  (b)  ooefBcieins 


BACKSCATTER/EXTINCTION  in  km ' 
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The  histognfiis  ahow  dut  die  nuahmiw  frequency  of  occiinence  occurs  at  an  extinction 
coefScient  of  about  0. 16  Imr*  and  at  a  badcscatter  coefSdent  of  about  0.01  km**. 


Fdr  inversion  of  the  vertical  lidar  measurementt  it  is  assumed  that  die  ratio  of  backscatter  and 


extinction  (P/a)  is  oMistant  diron^iout  the  mixed  layer.  Some  andmrs  even  apply  a  fixed  value,  hi 
our  opinion  diis  is  not  ccHiect  because  die  ratio  dqiends  on  die  aerosol  properties  whidi  vary  from 
day  to  day. 


Three  exanqiles  will  be  presented  of  periods  were  die  ratio  p/o: 
has  a  fixed  high  value 
has  a  fixed  low  value 
varies  in  time 


An  example  of  a  period  with  fixed  ratio  p/a  (but  varying  extinction  and  badcscatter  coefficients) 
is  shown  in  Figure  5.2.  In  die  period  from  10:00  to  12:(X)  ajot,  as  indicated  by  '1',  the 
measurements  were  performed  in  light  snow  resulting  in  enhanced  P  and  a  as  well  as  ahi^er  P/a 
ratio  than  in  the  period  thereafter.  In  the  period  after  12:(X)  ajn.,  as  indicated  by  '2',  there  was  no 
precipitation. 


TIME  of  tho  DAY 


EXTINCTION  in  km  ' 


Hgure  5.2:  Left  figure:  Backscatter,  extmctioa  and  ratio  p/a  on  January  15-16,  1985  (time  dwve  24 
bom  is  the  nest  day).  Rom  l(h00  a.iii..lo  12:00  am  it  was  snowing  UghOy* 
Right  figure:  scatter  diagrotiiowing  the  rdation  between  extinction  and  badcscatter. 


BACKSCATTER  in  km ' 
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The  in  Figure  show  that  after  12:00  a.m.  bodi  die  baducatter  and  die  extinction  varied 
by  more  than  a  factor  2,  but  the  ratio  varied  less  dian  10  %.  Fttrtfaennore  it  can  be  seen  diat 
die  ratio  P/a  in  snow  is  larger  dian  during  the  following  period  when  there  was  no  precqntadon. 
The  relation  between  P  and  a  has  also  been  visualized  in  a  scatter  diagnuD  (Hgure  5  Jd). 

Another  exan^le  of  a  period  with  a  iSxed  but  low  value  of  die  ratio  p/a  is  shown  in  Hgure  3.3. 
The  period  from  15:00  to  34:00  hours  (Le.  10:00  ajn.  the  next  day),  was  characterized  by  a  very 
good  visibility  (>40  km),  in  adiidi  p/a  drt^  from  0.22  to  less  than  0.01.  AMiough  the 
backscatter  and  die  extinction  coefScients  varied  more  dian  one  order  of  magnitude  during  die 
whole  day  (see  the  scatter  plots  in  Rguie  53b),  the  ratio  P/a  varied  about  20  %. 


EXTINCTION  in  km  ’ 


Figure  5.3:  Scatter  diagrain  of  badcscatter  versus  extinction  (a)  and  die  backscatter/extinction  ratio  as 
a  function  time  (b)  CO  December  11. 1985  (tinie>  24  hours  is  the  next  (taQr) 


Figure  5.4  shows  a  period  of  a  varying  and  slightly  decreasing  P/a  ratio  over  a  period  of  about  24 
hours.  Rrom  11:00  aJiL  to  07:00  pjn.,  die  ratio  p/a  varies  about  10  %  and  after  07:00  pjn.,  P/a 
decreases  gradually  from  about  0.9  to  about  0.6. 
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The  setected  results  show  that  the  ratio  fi/a  in  Figures  5.1,  5.3  and  S.4  are  couqparable.  The  low 
value  of  p/a  in  Hgure  5.2  differ*  from  die  odier  examples.  These  results  illustrate  diat  the  ratio 
P/a  is  not  a  fixed  value  but  may  vary  frmn  day  to  day  and  during  a  day. 


The  extinctum-  and  die  backscattei  coefficients  are  calculated  with  the  slope  mediod  from  the 
log-range  corrected  nignal  (linear  regressioo).  During  diis  process  also  the  standard  deviation  of 
the  CTrinction  and  of  die  backscatter  are  calculated.  The  standard  deviations  are  not  only  caused 
by  electronic  noise  but  also  atmosidieric  structures  (See  also  Kunz,  1992).  It  is  expected,  that 
due  to  the  atmospheric  structure  die  standard  deviation  decreases  with  die  length  of  the  selected 
interval  (more  sanqiles)  which,  in  turn,  dqpends  on  die  actual  extinction  (and  backscatter)  as 
explained  in  Appendix  A.  It  is  eiqiected  that  the  standard  deviation  reaches  a  minimum  for  a 
certain  value  of  extincdoo  (or  backscatter).  This  effect  has  been  observed  in  many  data  files. 
Hgure  S.6  shows  a  typical  example  of  this  effect  for  both  the  extinction  and  die  backscatter. 


EXTINCTION  m  Xw '  aACKSCATTEN  I"  Xm ' 


Figure  S.6:  Scatter  diagrams  of  the  standard  deviatioo  of  die  extinction  (a)  and  of  die  backscatter 
coefficienu(b)  as  a  Iboctioo  of  the  extinction  and  the  backscatter.  Scatter  diagram  (c)  shows 
die  relation  between  <$/p  and  doAx. 

The  analysis  of  many  of  the  data  files  has  dwwn  diat  die  relative  deviation  in  mtdiictum  (or 
backscatter)  has  indeed  a  V  shape  (See  also  Figure  3.1).  Furdieniioce  it  appears  diat  die  minima 
of  the  variation  in  both  die  extinction  and  die  badcscatter  occur  almost  simnltaiieooily. 
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relation  between  latincrinn  »firt  m 

In  the  Figures  52  to  5  J  aome  daily  periods  were  presented  of  tlie  ratio  P/a.  The  histograms  for 
the  extinction  and  the  hackscatter  are  presented  in  Figure  5.1.  It  is  also  interesting  to  review  all 
the  available  data  for  the  ratio  ^a.  However,  presenting  all  available  data  in  a  scatter  plot  of  P 
versus  a  yields  one  dark  qmt  Therefore  die  averaged  P  values  and  their  standard  deviations  were 
calculated  for  discrete  values  of  die  extinction  as  shown  in  Figure  S.7.  To  diis  end,  die  extinction 
has  been  partitioned  over  200  channels.  The  data  were  further  separated  into  periods  wdh  and 
without  prec^itadon.  Results  for  the  dry  situation  are  shown  in  Hgure  S.7b,  those  obtained  during 
rain  are  presented  in  Hgure  S.7c. 
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Hgure  5.7:  Reladon  between  extinctioQ-  and  backscatter  coeCBcients  for  all  data  records  (a),  for 
shuatkxu  without  piecqiitatico  (b)  and  with  precqiitation  (c).  The  \tppa  and  lower  traces 
correspond  with  the  mean  backscatter  plus  or  mimrs  one  standard  deviation. 

The  results  in  Hgure  5.7  indicate  a  non  linear  relation  between  P  and  a  ^^lich  becomes  less 
evident  for  exdnctims  larger  than  about  5  knr*.  This  is  probably  due  to  local  floctuatioiis  in  a 
and  Fbr  extinctunis  smaller  than  0.03  knr*,  inaccurate  determination  of  the  extinctioo  causes 
the  variaHlity  of  PAi,  see  Kunz  and  DeLeeuw,  1992.  Li  diat  case  the  backscatter  can  be 
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detennined  more  reliably  than  die  exdncdon.  Precq^ditiiHi  has  mly  been  lecwded  during  woddng 
hours  because  no  instrumoit  was  available  at  that  time  for  automatic  indication.  It  is  clear  from 
Figure  S.7  that  there  is  no  distinct  difference  in  the  overall  shape  of  the  three  figures.  This  means 
diat  (this)  lidar  provides  similar  results  bodi  in  dry  weadier  and  during  precqiitation.  For  low 
values  of  the  extinction  and  the  backscatter,  the  backscatter  seems  to  reach  a  limiting  value 
of  bout  0.004  km->  which  is  about  a  factor  3  larger  dian  die  Rayleigh  backscatter  value 
(1.3x10-3  knrO-  For  extinctions  larger  than  about  S  km‘^  a  second  mode  becomes  visible  in 
which  the  backscatter  is  about  one  order  of  magnitiiHf;  lower.  A  more  detailed  analysis  shows  that 
this  is  only  present  in  a  few  files;  thus  far,  die  reason  for  diis  behaviour  is  not  understood. 

relation  between  relative  humidity  and  extincritm  or  barWitgatter 
The  relation  between  the  relative  humidity  and  the  backscatter,  the  extinction  and  the  ratio  ^/a  is 
an  other  subject  of  interest  (e.g.  DeLeeuw  et  al,  1986  and  Tonna,  1991)  because  the  relative 
humidity  tends  to  increase  with  height  in  the  atmospheric  boundary  layer.  It  is  generally  assumed 
that  there  is  a  positive  correlation  between  the  relative  humidity  and  the  extinction,  the 
backscatter  and/or  die  backscatter/extinction  ratio.  From  a  review  of  aU  the  da^a  files,  it  i^ipears 
that  die  correlation  between  the  relative  humidity  and  the  backscatter  is  better  than  die  correlation 
between  the  relative  humidity  and  the  extinction.  The  relation  between  backscatter  and  relative 
humidiQr  has  in  many  cases  been  observed  to  loop  back  during  a  period  of  rising  and  falling 
relative  humidity  (cf.  Figure  S.IO).  lliis  may  be  caused  by  the  anomalous  growth  of  hygrosct^ic 
aerosol  which  has  beoi  reported  e.g.,  by  Winkler,  1971. 

Some  typical  examples  of  die  relation  betwem  the  lidar  measured  extinction-  and  backscatter 
coefBcients  and  the  relative  humidity,  are  shown  in  the  Figures  S.8  to  5.10.  Hgure  5.8  shows  the 
scatter  plots  of  the  extinction  and  the  backscatter  coefficients  versus  relative  humidity,  as 
recorded  over  a  24-hoiir  period.  These  plots  illustrate  diat  the  variation  in  the  extinction  at  RH>90 
%  is  larger  than  die  variation  in  the  backscatter.  The  backscatter  increases  dightly  when  the 
relative  humidity  increases  from  54  %  to  90  %.  For  higher  values  of  the  relative  humidity,  the 
backscatter  increases  much  stronger.  This  effect  has  been  obsoved  in  many  data  files.  However, 
we  noted  that  die  relaticm  between  backscatter  and  relative  humidity  varies  from  day  to  dtty. 
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Scatter  diagram  of  the  extinction  versus  relative  humidity  (a)  over  a  period  of  24  hours  on 
February  20, 1985  and  the  scattn  Vagram  of  die  backscatter  and  die  relative  humidity  (b). 
Note  that  the  variation  in  the  backscatter  is  smaller  than  die  variation  in  the  extinction. 

Figure  5.9  shows  another  scatter  diagram  of  die  backscatter  versus  the  relative  humidity,  recorded 
on  March  26,  1985  when  higher  values  of  the  relative  humidity  occurred.  The  measurenrent 
period  was  about  24  hours.  For  values  of  the  relative  humidity  smaller  than  97  %,  die  backscatter 
varied  by  a  Uctar  of  about  2,  but  for  higher  relative  humidities  the  variations  were  up  to  diree 
orders  of  magnitude.  It  is  noted  diat  accurate  relative  humidity  measurements  near  100  %  are 
extremely  difficult  Therefore,  results  from  routine  measurements  should  be  carefblly  interpreted 
in  this  regime.  When  the  recorded  relative  humidity  approaches  100  %,  diis  might  coire^KMid  to  a 
situation  where  the  relative  humidity  goes  slightly  over  100  %  (doe  to  inherait  instrument 
inaccuracy)  and  aerosols  can  be  activated.  The  activated  aerosols  grow  in  size  very  fast,  resulting 
in  a  strong  increase  in  both  die  extinctitm*  and  backscatter  coeffidents. 
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Figure  S.8: 


RELATIVE  HUMIDITY  In  % 

Figure  5.9:  Relatiao  between  beckscatter  and  relative  humidity,  for  high  values  of  the  RH,  over  a  period 
of  24  hours  on  Match  26,  198S.  For  RH<90%  diere  is  a  positive  correlatioD;  for  hitter 
values  of  the  RH  there  is  no  relatioa  at  all 


Figure  S.IO  shows  an  example  where  variations  in  the  relative  humidity  result  in  a  liysteresis 
loop'  in  the  backscatter.  This  bduviour  has  been  observed  many  times.O 
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Rguie  5.10:  Relation  between  the  backscatter  and  the  relative  humidity  over  a  period  of  24  hours  on 
February  7, 1985.  Note  the  cod  between  74  and  82  %  relative  hunudiM. 
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-  idation  between  visflalitv  and  «.rtiiictionw 

The  visibility,  die  backacatter  and  die  extinction  coefBcients  axe  an  evident  set  of  jMmxieters 
\^iich  shoukl  oonelate.  Results  of  this  kind  have  been  published  and  discuaed  in  the  past  (among 
odiers,  by  Lamberts,  1978,  Lamberts  and  DeLeeuw,  1986  and  Bertolotti,  1969,  1978).  Some  of 
the  relations  derived  from  our  data  set  will  be  presented  below. 

hi  Hgure  S.ll  an  exaiqile  is  shown  of  a  24  hour  period  triieie  bodi  die  extinction  and  the 
backscatter  coefBcients  are  well  conelated  with  the  visibility.  For  high  visibilities,  die  deviatitm 
in  the  extinction  is  latger  dian  the  deviation  in  die  backscatter. 


VISIBILITY  in  km 


Figure  S.ll:  Scatter  diagram  of  die  extinction  coefficient  (a)  and  of  the  backscatter  coefficient  (b)  as  a 
function  of  the  visibility  over  a  period  of  24  hours  on  March  10,  lSt86  (time  above  24  hours 
is  the  next  day).  Note  dist  both  die  extinction  and  the  backscatter  correlate  well  widi  die 
visibility. 


Anatysis  of  all  validated  data  files  diows  diat  die  relation  between  visibility  and  backscatter  or 
between  visibility  and  extinction  is  not  always  as  good  as  shown  in  Hgute  S.ll.  hi  many  cases, 
the  visibility  anrelates  either  better  widi  die  backscatter  or  with  the  extinction,  hi  odier  cases  no 
correlation  at  all  is  found.  Figure  5.12  shows  an  exan^ile  where  die  visibility  conelates  b^ter 
with  the  extinction  (range  interval  fiom  about  0.5-2  knr*)  than  with  die  backscatter.  hi  Figure 
5.13,  however,  where  die  extinction  varies  ftocr  about  0.03-1.5  km*>,  die  visibility  correlates 
better  with  the  backscatter. 
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Hgure  5.12:  Scatter  diagmn  of  die  extinction  (a)  and  the  beckacang  (b)  as  a  function  of  the  visftrility 
measuied  ovg  a  24  hour  period  on  April  7,  1986.  Note  that  die  conelation  between 
extinction  and  visibility  is  better  dum  the  coneladoo  between  backscatter  and  visibility. 
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Kguie  5.13:  Scatter  diagram  of  the  extinction  (a)  and  the  backacalter  (b)  as  a  function  oi  the  visibility 
measured  ovg  a  24  hour  period  on  ^ril  18.  1986.  Note  diat  die  conelation  between  die 
backscaltg  and  the  visibiliQ'  it  better  dian  the  conelmion  between  die  extinction  and  the 
visibility. 


The  last  exanqile  in  this  series  shows  diat  dien  are  also  periods  where  neidier  die  backscatter  nor 
the  extinction  correlate  well  widi  die  visibility.  This  is  shown  in  Rgure  5.14. 
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Figure  S.14:  Scatter  diagram  of  the  extmctioD  (a)  and  the  backacattcr  (b)  at  a  finctioo  of  Ae  visibility 
measured  over  a  24  hour  period  oo  FObiuaiy  21,  1986.  Exan^  of  a  poor  coirelatioo 
between  die  visiMlity  and  die  scattering  coefficients. 

S.2  Analysis  of  vertical  measurements 

Knowledge  of  dm  vertical  stmctuie  of  die  atmosidiere  is  of  itE^portance  for  many  disciplines. 
Especially  processes  in  die  atmosidieric  boundary  layer  (iq>  to  1-2  km)  ^y  an  imponant  role  for 
human  life.  These  indode  both  observable  processes,  like  haze  and  wind,  and  the  transport  of 
pollution,  mixing  of  air  from  difSerent  altitudes,  etc.  The  vertical  or  slant  path  visibility  is  an 
inqiortant  parametCT  in  avionics;  the  altitude  and  behaviour  of  the  mixed  layer  ate  also  important 
for  meteorology. 

Results  from  measmements  of  vertical  atmoi;dieric  extinction  profiles  during  die  VISA  project 
are  presented  for  selected  periods.  As  far  as  possible,  exanqiles  are  presented  of  profiles  vihich 
were  measured  during  several  di^s  to  one  week  in  each  mondi  of  the  year,  to  give  an  impression 
of  possible  seasonal  effects.  Results  ate  presented  as  fidse  cdor  i^its,  v^iere  die  altitude  is  plotted 
along  the  vertical  axis  and  the  time  is  plotted  along  die  horizontal  axis.  The  extinction,  vriiich  is 
the  dihd  parameter  in  these  figures,  is  coded  in  folse  cdor.  The  color  scale  is  plotted  <m  top  of 
each  figure  for  reference.  Hourfy  intervals  are  indicated  by  diin  vertical  tics  at  die  top  and  die 
bottom  while  die  days  are  separated  by  white  vertical  twnm. 
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-  rBBiltoffimMnrl9!t^ 

Vertical  extinctioo  profika,  as  measorad  in  the  period  JaDuaiy  21-24,  1985,  are  presented  in 
Figoe  5.15.  An  external  minor  was  used  to  leadi  a  maximum  altitude  of  2000  m.  Hie 


January  1985 

Figure  5.15:  Vertical  extinctioa  profiles  coded  in  fiUse  color  from  January  21  (8  a.ni.)  to  January  24 
(8  a.in.),  1985.  Hie  maximom  altitude  was  2000  m. 


On  Mondqr  January  21, 1985,  from  08:21  ajn.  to  about  11:30  am,  the  extinctioo  profiles  could 
only  be  measured  to  an  altitude  of  about  80  m  due  to  a  ground  fog.  Hie  extinction  coefficient  was 
very  high  (about  2  knr^).  Between  11:K)  am  and  04:00  pm,  die  extinction  at  ground  level 
decreased  and  prtdiles  could  be  measured  up  to  an  altitude  of  about  14(X)  m  Hien  a  cloud  layer 
descended  from  about  500  m  to  about  200  m  at  07:(X)  pm  Subsequent^  die  ckud  layer  rose 
gradually  to  diout  300  m  and  die  extinction  at  ground  level  decreased  cooespondingly  to  about 
0.2  km*!.  After  muinighf,  •  new  clood  Imyer  was  detected  at  about  1200  m  A  low-aldtnde  cloud 
base,  decreasing  from  about  300  m  altitude  to  about  200  m  altitude,  was  present  until  about  ll.-OO 
am.  (Ml  January  22.  At  diat  mnnent,  a  doud  layer  was  detected  at  about  500  m  altitude  during 
about  two  hours;  die  wind  veloci^  decreased  from  about  8  m^  to  about  4  m/s.  Hiis  doud  layer 
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cuddenly  dia^ipeaied  after  01:00  pm  but  odwr  cloud  leflecrions  were  obaerved  at  altitudes 
betweea  about  1400  m  and  1600  m.  Periotto  with  clear  aides  occmed  during  die  nigbt  (human 
observation),  accompanied  by  relatively  low  extinction  values  at  ground  level  After  mid-nig^t, 
the  extinction  at  ground  level  incieaaed  until  about  05M)  am  on  January  23.  Between  09:00  am 
and  Q2K)0  pm  die  cloud  base  varied  between  500  and  1500  m  altitude.  Starting  aroui^  noon  diis 
day,  patdies  of  douds  were  delected  at  different  altitudes  during  the  afternoon,  evening  and  ni^t 
The  wind  was  strong.  On  January  24  the  bducvionr  was  similar  to  that  observed  during  the  night 
before. 


PiitgniMitwi: 

The  results  from  Figure  5.15  show  die  dissolution  of  die  ground  haze  in  the  morning  of  January 
21  and  the  increase  in  die  depdi  of  die  mixed  layer.  The  altitude  of  die  cloud  layers,  whidi  are 
indicated  1^  small  white  vertical  tines  or  small  white  patches,  was  very  instable  and  varied 
between  80  m  and  2(XX)  m  during  dus  period.  Note  that  the  inSnence  of  clouds  below  S(X)  m 
altitude  is  percqitible  in  the  extinctkm  some  hundred  meters  below  die  cloud  base. 

*  example  of  Febniatv  1985 

Figure  5.16  shows  extinction  profiles  to  an  altitude  of  2000  m  over  a  period  of  5  consecutive 
days.  The  measurements  started  on  February  4,  at  08:43  am  and  ended  on  February  8,  at  03:28 
pm 


February  1985 

Figure  S.  16:  Veitical  aamctioii  profile*  coded  in  Abe  color  in  a  height  versus  time  plot.  Fniod:  February 
4, 08:43  a.m.  to  February  8, 1985, 03:28  pjn. 


At  die  beginning  of  diis  period,  some  gioiind  haze  was  present  In  the  log  was  recorded  dial  the 
sky  was  covered  with  thin  cirrus  (not  visible  in  die  figure).  Starting  at  11:00  am.,  a  layer  with 
increased  extinction  was  detected  at  750  m  altitude.  This  was  probably  not  a  cloud  layer  because 
there  was  also  sunshine  during  diis  period.  The  extinction  at  ground  level  decreased.  At  05:00 
p.iiL,  the  depdi  of  die  mixed  layer  started  to  decrease  from  about  750  m  altitude  to  about  200  m  at 
08:(X)  am.  on  February  5.  (This  resulted  in  a  ground  haze,  die  extinctiMi  at  ground  level  increased 
until  about  07KX)pm.)  This  ground  haze  sti^ed  during  die  night  and  moming  (observed  visually). 
During  die  ni^  die  sky  was  clear  and  diere  was  some  frori.  After  08:00  am.  on  Frivuary  5,  die 
mixed  layer  rose  gradually  to  aboid  7(X)  m  altitude,  udiich  was  reached  around  10:(X)  pm  The 
extinction  increased  and  douds  were  visually  observed  (not  shown  in  the  figure).  Aldiough  the 
altitude  of  die  mixed  layer  was  constam  during  the  night  and  the  moming  on  February  6,  die 
of  die  mixed  li^  diowed  a  strong  coovdving  bdwviour  widi  a  perkid  of  almost  two  hcnirs.  The 
extinction  at  ground  level  increased  significantiy  resulting  in  a  misty  moming.  Between  around 
07K)0  am.  and  08K)0  am,  die  height  of  die  fog  layer  extended  to  about  100  m  altitude.  I¥om 
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08:00  ajn.  to  03:00  pjn.  oo  Febniaiy  6,  a  cloud  layer  was  detected  at  about  1500  m  altitude. 
Around  01:00  pro.,  die  fog  li^  qipeated  to  $epmtic  into  a  ground  fog  layer  and  a  clood  rising 
from  about  400  m  to  about  1000  m  altitude.  Rtom  06:00  pro.  to  about  lOKX)  pro.,  a  cloud  layer  is 
detected  ^riiich  rose  fimn  about  100  m  to  about  1300  m  altitude.  During  tbe  nigbt,  tbe  dqpdi  of  the 
mixed  h^er  extended  to  about  1300  m  altitude  while  a  slight  variatioo  in  extinctioo  was  visible; 
ten^eratnre  decreased  below  the  fieezing  point  On  Fdxuary  7,  bttween  06:00  and  07:00  aro., 
die  height  of  die  mixed  layer  decreased  rapidly  from  1000  m  to  100  m.  After  08:00  aro.  the 
mixed  layer  rose  gradually  from  about  100  m  to  about  800  m  at  09:00  pro  During  die  night  and 
the  morning  oa  Fdiruary  8,  die  height  of  die  mixed  layer  varied  somewhat  around  500  m  altitude 
^iriiile  the  extinction  itself  decreased  gradually  to  a  minimum  around  07:00  aro.  After  00:30  pro, 
a  npdly  descoiding  cloud  layer  was  detected. 


BMCBSaOB 

The  moming  warming  effect  is  clearly  visible  on  Fdvuary  5, 7  and  8  starting  at  about  09:00  aro. 
The  effect  of  increaang  extinction  during  die  night  is  visible  on  February  4  and  5  and  partly  on 
February  7.  The  depdi  of  die  mixed  layer  varied  from  less  dian  100  m  to  more  dian  2  km  and  is 
the  average  dqidi  was  about  700  ro  Tbe  meandering  effect  of  die  top  of  die  mixed  layer  is  clearly 
visible  in  the  night  fimn  Februsy  5  to  6  and  in  the  moming  of  February  6. 

-  example  of  March  1985 

Hgure  5.17  shows  an  exanqile  of  extinctitm  profiles,  coded  in  false  color,  as  measured  with  lidar 
during  five  consecutive  days  fixmi  March  4, 08:45  aro,  to  March  8, 03:30  pro 
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^tarch  1985 

Figure  3.17:  Vertical  extinctkai  pn^es,  coded  in  false  color,  from  March  4-8.  1983.  The  period  from 
8:30  ajD.  to  02:30  pjn.  on  March  4,  wu  diaractetized  by  rain  showers.  The  strong 
extincdon  during  Miu^  6  and  7  was  caused  by  a  ground  haze. 

The  period  from  08:30  aon.  to  02:30  pjn.,  c»  March  4,  was  dominated  by  many  rain  showers. 
The  cloud  base  was  detected  around  200  m.  After  02:30  pm,  die  cloud  base  rose  to  about  900  m. 
After  04:30  am,  die  clouds  vanished  and  the  son  ai^eared.  During  the  clear  night,  the 
tenqieratore  at  ground  level  decreased  below  0*  C  and  ground  haze  was  formed.  The  height  of  die 
mixed  layer  remained  below  250  m  altitude.  The  extinctkm  started  to  increase  at  about  10:00  pm 
on  March  4  until  08:00  am  on  Manh  5.  At  about  08:00  am,  die  ground  haze  dissolved  and 
disqipeaied  con^letely  at  01:00  pm  From  diat  time  on,  die  extinctkm  at  ground  level  inoeased 
strongfy  and  reached  its  maximum  value  at  midnigjit  The  depth  of  die  mixed  h^er  remained 
almost  ccmstant  at  an  altitude  of  about  100  m.  At  01:00  pm  cm  March  6,  die  ground  haze 
distfipeared  and  the  sun  was  visible  during  about  twp.hours.  This  was  followed  by  an  earty 
morning  ground  haze  which  lasted  until  about  11:30  am  on  March  7.  (The  high  extinction  vahies 
are  not  real  but  are  causes  by  singularities  in  the  forward  inverskm  method.)  During  die  ni^t 
from  March  6  to  7,  die  dcy  was  clear.  On  Mardi  7  at  11:30  a.m,  die  mixed  layer  deepened  until 
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about  04:00  pjn.  From  Otat  time  to  about  09:00  a.m.  on  Mardi  8,  die  dqidi  of  die  mixed  layer 
meandered  over  an  interval  of  about  SO  m  around  an  altitude  of  about  100  ol  At  09:00  a  m.  a 
cloud  was  detected  and  short  dureafier  a  new  ground  haze  fonned. 

PaCTBMQO 

Outing  the  March  period,  diere  was  naich  ground  fog  widi  a  thickness  of  about  100  m.  Diumal 
effects  were  visible  in  the  night  from  March  4  to  S  and  during  die  moming  of  Match  7.  Several 
tunes,  the  forward  inversion  method  failed  in  die  ground  fog  (Match  5  from  10:00  pjn.  to  March 
7  at  11:00  ajn.)-  During  modest  fog  (March  7,  06:00  to  08:00  ajn.)  die  vertical  atmo^heric 
extinction  profiles  could  be  measured  and  inverted.  This  can  be  seen  from  die  structure  at  larger 
altitudes.  On  March  7,  the  mixed  layer  rose  with  a  velocity  of  about  SO  m  per  hour. 

example  of  April  results 

The  results  of  the  vertical  extinction  profiles  measured  from  April  IS  to  April  19,  198S,  are 
shown  in  Hgure  S.18.  This  dry  period  was  characterized  by  a  clear  sky  and  ground  haze.  On  April 
IS  at  08:S0  am.,  the  measurements  started  under  a  cloudy  sky  with  a  cloud  base  at  about  500  m. 
After  09:00  am.,  die  extinction  in  die  mixed  layer  decreased  until  about  1 1:00  pm.  At  01:(X)  am. 
tm  April  16,  a  cloud  was  detected  at  about  2S0  m  and  the  extinction  at  ground  level  increased 
rapidly.  At  06:00  am.,  a  new  cloud  layer  was  detected  at  about  7(X)  m  altitude  which  remained 
there  until  about  08:00  pm.  Note  that  die  extinctioa  below  die  cloud  started  to  increase  already  at 
an  altitude  of  about  3(X)  nu  Betwem  05:00  pm.  and  07KX)  pm.,  a  haze  layer  was  observed  to 
move  in  from  the  sea.  This  resulted  in  an  increase  in  die  extinction  between  200  and  600  m 
altitude.  Doting  die  night  from  .^pril  16  to  .^itil  17,  a  ground  haze  develt^ied.  The  depdi  of  this 
layer  decreased  slowly.  Minimum  altitude  was  reached  between  03:(X)  and  04:00  am.  Rom  that 
dme  on,  die  top  of  the  ground  fog  layer  ascended  from  about  100  m  to  about  250  m.  A  second 
li^er  was  detected  vdiich  rose  fiom  abmit  200  m  to  SSO  m  aldtode.  Rom  01:00  to  03:00  pm.,  die 
atmoqihere  above  100  m  was  clear.  At  about  04:00  pm.,  cumulus  clouds  were  observed  and  at 
0S:(X)  pm.  a  new  ground  fog  layer  developed  which  stayed  until  01:00  pm.  on  April  18.  At  dut 
tiine  die  extinction  decreased  due  to  deqieaing  of  die  mixed  layer.  From  about  04:00  pm.  on 
April  18  to  about  02:(X)  am.  on  ^iril  19  the  diickness  of  the  mixed  layer  decreased.  Subsequently 
a  new  layer  was  detected  at  an  altitude  of  about  700  m  which  descended  slowty  until  10:00  am. 
to  an  altitude  of  about  KKlntFrcnn  that  time  cm  dds  layer  ascended  and  dissolved. 
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Figure  S.18:  Vatkal  extnictioD  profiles,  coded  in  fidse  color,  in  die  period  April  IS  to  19,  198S; 

maximum  altinide  1000  m.  The  mixed  layer  varied  in  ahhude  from  about  400  m  to  more 
than  1000  m  altitude.  No  precipitation  was  observed  during  this  period.  Oround  was 
present  fiom  April  Id,  02:00  a.in.  to  ^xil  18, 12:00  ajn. 


Discussion: 


Hie  Ajoil  period  showed  a  mixed  layer  with  a  varying  thickness  fiom  about  200  m  to  more  tfian 
1000  m.  The  lypearance  and  dissohidoa  of  fog  layers  were  clearly  mapped.  In  part  of  the  penod, 
die  inversion  of  some  lidar  returns  failed  due  to  die  very  strcmg  extinction  at  ground  leveL  During 
two  days,  a  rise  of  die  mixed  li^  in  die  moming  has  been  observed.  The  rate  was  about  120  m 
per  hour.  The  data  fiom  figure  5.18  also  show  diat  die  extinction  undenseadi  a  cloud  can  extend 
over  several  hundreds  of  meters. 


example  in  May  1985 

Results  of  vertical  atmospberic  extinction  profiles  as  measured  with  lidar  in  dteprec^ntatkm-fiee 
period  on  May  6-10, 1985  are  shown  in  figure  5.19.  Hie  measurements  started  on  May  6, 1985  at 
09:23  ajn.  The  log  reported  diat  cinus  douds  covered  die  dcy  and  diat  some  ground  haze  was 
present  wliidi  stayed  until  about  03:00  pjn.  The  thickness  of  die  mixed  Itya  grew  while  die 
extinction  at  die  tarface  remained  unchanged.  Also  udien  die  son  came  dnongh,  die  surface 
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extiiicti(Hi  did  ncrt  decrease.  During  die  evening  and  the  night  (including  die  morning  of  May  7) 
the  ground  extinction  increased  and  reached  maxima  armmd  02:(X)  ajn.  and  08KX)  ajn..  At  08:00 
ajn.,  the  sky  was  fully  overcast  (which  is  not  visible  in  the  lidar  results).  From  09:00  ajn.  to 
03:00  pjn.  the  height  of  the  mixed  b^er  decreased  somea^iat  and  after  03:00  pjn.  it  increased 
again.  The  extinction  at  ground  level  also  increased  and  readied  a  maximum  at  about  10:00  pjn. 
Subsequently  the  ground  exdnctioo  varied  somewhat  and  reached  a  new  maximum  at  06:00  ajiL 
(m  May  8.  At  about  07:30  ajn.  die  altitude  of  die  mixed  layer  started  to  rise,  reaching  its 
maximum  at  about  04:00  pjn.,  udiile  die  ground  extinction  decreased.  Than  a  cloud  layer  was 
detected  at  about  1000  m  altitude,  which  disappeared  again  at  08:00  pjn.  From  that  time  die 
ground  haze,  widi  a  thickness  of  about  100  m,  was  detected  and  lasted  for  the  rest  of  die  period. 
From  about  10:(X)  pjn.  on  May  8  to  09:(X)  ajn.  on  May  9,  the  surface  haze  layer  was  ciqiped  by  a 
clear  layer.  On  May  9  at  09:00  am.  the  sky  was  clear  with  sun  shine.  After  09:00  am.  a  new 
ground  haze  developed  and  this  layer  rose  widi  a  velocity  of  about  100  m  per  hour.  The  fog  was 
also  present  during  M^r  10  which  resulted  in  a  limited  range  of  die  lidar. 


May  1985 


Rgure  5.19:  Veitical  extinctiaa  fnofiles,  coded  in  fidse  cdour,  in  a  time  versus  hei^  (0-1000  m)  plot 
Period  Ctom  May  6  to  May  10, 1985.  The  period  is  characterized  with  intervals  of  mist  and 
haze  but  without  pracipitatiaii. 
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Discussion: 

Dtumal  variations  are  visible  on  May  6,  7,  8  and  9.  On  May  7,  8  and  9,  the  development  of 
ground  haze  has  been  observed  widi  different  extinction  values. 


result  of  Julv/Augnst  1985 

The  vertical  atmosiAeric  extinctioa  profiles  measured  during  die  period  July  29  to  August  2, 
1985,  are  presented  in  figure  5.20.  During  July  31  the  sky  was  fiiUy  overcast  without 
precipitation.  The  height  of  the  cloud  base  varied  from  500  m  to  over  1000  m  but  die  cloud  layer 
could  not  be  distinguished  clearty.  During  the  night  on  July  29  and  31,  ground  mist  layers  were 
detected  which  reached  altitudes  of  respectively  250  and  500  m.  A  diurnal  effect  is  visible  in  die 
morning  and  in  the  evening  of  July  31. 


July  1985  /  August  1985 


Figure  5.20:  Vertical  extinction  profiles  coded  in  frlse  color,  in  a  height  versus  time  plot, 
measured  with  lidar  in  the  period  July  29  to  August  2, 1986. 
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CTMipte  gf  OctebCT  1985 

llie  vertical  atmoqphedc  extioctum  profiles  measured  daring  the  fourdi  wedc  in  October  1985, 
are  presented  in  Hgore  5^1.  No  clouds  were  detected  in  diis  period.  During  die  first  diree  days, 
die  height  of  die  mixed  layer  varied  between  about  300  m  and  about  600  m.  b  die  last  two  days, 
the  depdi  of  die  mixed  layer  was  about  2SO-3SO  iil  During  die  night  from  21  to  22  October  and  in 
die  ea^  morning  of  October  22,  the  dqith  of  the  mixed  layer  decreased  to  about  300  m  and  die 
extinction  at  ground  level  increased.  Around  09:00  ajn.,  die  height  oi  the  mixed  layer  rose  again 
to  about  600  m  with  a  velocity  of  about  40  m/hour  and  die  extinction  at  ground  level  decreased. 
At  about  07:00  pjn.,  die  mixed  layer  remained  at  a  constant  level  of  about  600  m  until  at  least 
midnight  on  October  23.  (Due  to  a  system  error,  die  lidar  was  not  in  operation  from  midnight 
October  23  until  about  09:00  ajn.  diat  day.)  After  06:00  pjn.  on  October  23,  die  height  of  die 
mixed  layer  decreased  significantly  with  a  rate  of  about  80  m/hour.  The  extinction  at  ground  level 
also  decreased  and  reached  a  minimum  value  at  midnight  Note  that  this  is  opposite  to  the 
behaviour  during  the  night  from  Octoba  21  to  22.  During  the  morning  of  October  24,  die 
extinction  at  low  altimdes  increased  imril  about  08:(X)  ajn.  The  mtoimum  value  of  die  ntinction 
occnnred  at  about  03:00  pjn.  and  increased  gradually  until  October  25,  09:00  a.ni.  At  diat  time, 
die  height  of  die  mixed  layer  started  to  rise. 
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Figure  S.21:  Vertical  extinctioa  profiles,  coded  in  false  colour,  in  a  time  versus  height  diagram.  Period 
October  21  to  October  25, 1985.  Varying  mixed  layers,  both  in  altitude  and  in  extinction  can 
be  distinguished. 


PiggwgiQB; 

The  vertical  extinction  profiles  measured  in  die  October  period  clearly  show  the  process  of  a 
developing  mixed  layer  widi  increasing  extinction  at  low  altitudes  from  midnight  until  about 
09:00  am. 

-  evamnle  nf  Nnvemhta-  Hate 

A  soies  of  vertical  atmospheric  extinction  measurements  in  November  1985  is  presented  in 
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November  1985 

Figure  5j22:  Votkal  exdnctiaii  profiles,  coded  in  fidse  cokw,  in  a  time  versus  hei^  diagnm.  Period 
November  4  to  November  8, 1985. 

On  November  4,  die  log  reported  diat  the  sky  was  fully  overcast,  but  the  cloud  reflections  are  not 
visible  in  Figure  5.22.  (Le.  Cloud  base  was  hi^er  dian  1000  m.)  During  die  early  morning  on 
November  5  the  cloud  layer  descended  n^iidly  with  a  rate  of  about  ISO  m/hour  and  a  fog  layer 
with  a  diickness  of  about  200  m  developed  whidi  remained  until  about  07:00  ajn.  Dhitnal  effects 
ate  slightly  visible  in  die  extinction  at  low  altitudes  in  the  evening  of  November  4  and  in  the 
morning  of  November  7  and  8. 

Ptscussion: 

This  Novendier  period  is  mailced  by  relatively  low  extinction  values  at  ground  level  (excqit  in  die 
morning  of  S  November)  and  a  variable  extinetkm  both  in  altitude  and  in  time.  Qouds  were 
detected  at  different  altitudes  at  irregular  moments.  The  mixed  layer  varied  between  about  250  m 
and  1000  m  and  more. 


fitaMwl  rfi«CMiii>inn 


False  cokr  piescntarion  of  die  vertical  exdnctioD  pn^iles  is  a  suitable  m^bod  to  peeseut  large 
amounts  of  qidar)  data  in  a  quasi  daeediincasiooal  figure.  Although  die  figures  give  a  qualitative 
inqxessuMi  of  die  tenqwnd  and  spatial  behavioar,  die  extinction  values  are  quantitatively 
available.  By  reviewing  the  data  files,  it  is  noted  that  in  contrast  to  what  is  generally  assumed,  die 
increase  in  die  dqidi  of  the  mixed  layer  during  die  moming  is  not  alwi^s  dearly  visibk  in  the 
measured  data.  This  is  also  die  case  for  the  descend  of  die  mixed  layer  in  the  evening  and  night 
The  reason  for  diis  might  be  the  fint  dial  many  fronts  with  varying  wind  direction  and  different 
air  mass  (with  different  a  and  P)  the  atmoqilietic  conditions  in  die  Nedierlands  (See 

e.g.  Lamberts  and  OeLeeuw,  1986).  Fdrdieimare,  turbulent  mudng  due  to  solar  heating  was  in 
many  cases  reduced  by  die  presence  of  clouds.  If  presmt  die  ascend  rate  of  the  top  of  die 
boundary  layer  is  on  the  order  of  1(X)  meter  per  hour,  as  derived  fiom  die  figures  presented. 
Variations  in  altitude  of  the  mixed  have  clearly  been  nuqiped  widi  the  lidar.  The  hei^t  of 
die  mixed  b^er  varies  generally  slowly  (over  hours),  but  also  sudden  changes  on  tiine  scales  of 
scxne  minutes  (including  abnormal  conditirms  like  rain)  have  been  observed.  From  die  available 
measurements,  it  is  possdile  to  detennine  the  padi  integrated  extinction  bodi  in  vertical  smise  as 
well  in  slant  padi. 

S.3  Statistical  review  of  die  vertical  extinction 

The  frequency  of  occurrence  and  the  cumulative  distributions  of  the  extinction  at  altitudes  of  IS, 
140, 265,  390,  S2S,  640, 76S,  890  and  lOlS  m  above  ground  level  have  been  determined.  This  is 
at  intervals  of  12S  m  phis  the  altitude  of  die  system  at  15  m.  To  sk^  anomalies,  mly  profiles  in 
periods  without  precqiitation  have  been  incloded  in  diis  analysis. 

-  wmMncdrewte  frwa  and  L9S^ 

The  histograms  of  die  extinctioo  at  the  above  mentioned  altitudes,  derived  fimn  all  die  1985  and 
1986  data  files,  are  showr  in  Figure  5.23.  The  cumulative  distributions  are  shown  in  Figure  5.24. 
The  mean  vafaie  and  the  standard  deviation,  both  derived  from  the  cumulative  distribotioo  and 
plotted  as  a  frtnction  of  height,  are  presented  in  Rgure  5.25.  The  value  of  0.2  km-^  at  15  m 
altitude,  is  equal  to  the  results  of  Figure  5.1  (frequency  distiibotiim  of  die  horiziHitBl  exdnctiao). 
The  high  frequency  of  oocuttenoe  in  die  fiigure  at  IS  m  and  140  m  me  due  to  die  dioice  of  a 
bonndaiy  condition  in  diose  sitnadmis  were  die  horizontal  measurement  did  not  provide  an 
adequate  solutuin. 
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Tbe  biatograms  in  Rguie  523  have  diffeiem  duqwa  at  altitudes  below  and  above  400  m.  Below 
dus  abitude,  a  itmltmiode  bduvioor  is  diaerved;  above  diis  abitnde,  die  frequency  distributkais 
qiproacb  a  skew  Poisson  distnbatioii  widi  a  maTwruw  Bkelihood  at  about  0.02  km**.  Tbe  widths 
of  tbe  distributiaos  see  moce  than  1 S  decades,  Veiy  large  extinctions  (above  1  knrO  are  rare.  The 
histograms  at  390  m  and  higher  shift  gradual^  to  someiAat  lower  extinctkm  values.  This  means 
diat  dieie  is  a  slight  tendency  diat  die  extinction  decreases  as  the  altitude  increases.  This  dSect  is 
better  visitde  in  Rgure  5.25. 
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Figure  S.2S:  Mean  (solid  squares)  and  standard  deviadoo  (open  circles)  oi  tbe  averaged  estmcdon 
as  a  fimctkn  of  altitude  from  all  1985  and  1^  data  files  (11,000  profiles)  widiout 
rain. 


Id  1986,  a  more  accurate  predptation  registration  was  established  and  die  wiy  of  data  processing 
was  modified  somenbat  Therefore  die  1986  results  (5(X10  profiles)  were  analyzed  sqwraiely  (for 
periods  widiout  precqritation).  The  frequency  and  die  cumulative  distributions  are  shown  in 
Rgnies  526  and  527.  In  Rgure  5.28,  die  mean  and  die  standard  deviation  are  presented  as  a 
functkn  of  height 


ngme  5^:  Diftribaticm  cX  die  atmoqtiiefk  eaffincrioo  at  idne  diffenni  alritmte*  firom  the  1986  data 

(3000  profiles).  Only  perimb  vddwot  pfec4>><><^  X* 
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FigmeS^:  CainulMiw  diWittiilioii  of  die  1986  data  m  nine  different  ihinidei.  See  «l»o  figure  5^ 
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The  fre<pieiicy  distribotioiu  in  Hgure  5Ji6  for  die  tow  altitudes  (<400  m)  are  diffoeat  from  diose 
in  Hgwe  SJ23.  Multi-modal  distribntkms  are  not  observed,  only  single  mode  distribtttkms.  A 
small  qnke  around  10-3  knr*  coneymds  to  the  default  extincticm  in  situattoos  (of  very  good 
visibility)  v^ien  die  lidar  provides  a  negative  extinction.  The  peak  in  the  frequency  distribution  at 
15  m  altitude,  around  0.2  km'>  coneaponds  widi  die  mean  extinction  as  shown  in  Hgure  5.1.  The 
mean  extinction  shifts  slowly  from  abtuit  0.2  knr*  at  20  m  altitude  to  about  0.08  km-i  at  one 
kilometer  altitude.  At  higher  altitudes,  die  extinction  is  qvead  over  a  larger  interval  which 
extends  mostly  towards  smaller  extinction  values.  The  averaged  values  and  die  standard 
deviations  of  the  1986  data  show  a  somewhat  rimiiar  hehaviour  as  die  combined  results. 
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Figure  5.28:  Mean  (solid  squares)  and  standard  deviatioo  (open  circles)  of  the  extinction  (1986  data;  5000 
proftles)  as  a  function  of  altitude. 
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6  DISCUSSIC»>IANDCX)NCLUSIONS 

During  a  period  of  19  mondu,  from  January  198S  until  August  1986,  vertical  atmo^Aeric 
extinctum  profiles  were  measured  widi  a  small  lidar  system  (critical  radar)  to  an  attitude  of 
maximum  2000  m.  The  requited  hqnt  boundary  condition,  for  the  (fraward)  inverskm  of  the 
signals,  was  measured  widi  die  same  lidar  system.  Data  has  been  collected  continuously  over 
periods  of  maximum  5  woridng  days.  After  validation  and  combinatiim  widi  die  actual  meteo 
data,  a  data  base  has  been  established  of  22,800  records. 

The  lidar  system  and  the  transfer  functions  of  die  electronic  conqxuients  were  absolutely 
calibrated.  Thus  both  die  badcacatter*  and  die  extinction  coefficients,  as  well  as  the  extinction 
profiles,  were  measured  quantitatively.  It  has  been  shown  diat  the  accuracy  of  the  trigger  moment 
of  the  recording  system  and  die  accuracy  of  die  transfer  functim  of  die  logaridimic  anqilifier  play 
an  inqiortant  role,  especially  in  situations  of  very  low  extinction  coefficients.  The  inverskm 
method  qiplied  has  been  validated  during  a  field  experiment  in  1983,  where  aerosol  and 
nefdieloiiieter  profiles  were  measured  simultaneously  with  die  lidar  (Lindberg  et  al,  1987). 

The  measured  «itinction  profiles  are  the  basic  information  for  nalmUring  the  vertical  or  slant 
path  visibility.  This  parameter  is  very  inqxntant,  e.g.,  for  meteorolti^,  aviation  and  for  die 
prediction  of  the  performance  of  electro-optic  sensor  systems.  At  this  moment,  diere  are  no  odier 
efficient  techniques  to  determine  the  vertical  optical  structure  of  the  atmoqihere. 

Analysis  of  the  data  base  has  led  to  die  following  crmclnsions: 

tt  is  estimated  diat  die  forward  inversion  method  provides  reliable  information  on  the 
atmoqdieric  vertical  extinction  jaofile  in  about  95%  of  die  cases.  Only  during  periods  of 
mist  or  ground  haze,  the  reverse  inversion  medwd  mi^  provide  better  results. 

In  general,  there  is  no  fixed  relation  between  die  extinction  and  die  baclczcatter  coefikients. 
Neverdidess  periods  can  be  distingniahed  where  a  constant  relaticm  is  clearty  present  despite 
the  variatum  in  die  badcacatter  and  the  extinction. 
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For  relittive  hmniditiea  aaudkr  dum  90  %,  both  tbe  badcacatter  aod  Ote  extinctk>n 
co^ficknU  increaae  somewiiat  with  the  relative  humidity.  For  larger  values  of  the  relative 
humidity,  die  incremeot  is  stronger  but  not  well  defined.  Due  to  the  struig  variabilhy  in  die 
measured  quantities,  which  mty^  be  caused  by  different  air  masses,  it  is  not  possible  to 
(ferive  an  empirical  relatitn  between  relative  humidity  and  die  extinctuMi  or  badcscatter 
coefBcknL 

hi  many  cases,  a  liy  steresis'  loop  has  been  observed  in  the  relation  between  tbe  badcscatter 
coefficient  and  die  relative  humidity.  We  presume  diat  diis  is  caused  by  die  shift  in  tht 
partide  size  distribotions  in  response  to  changes  in  relative  humidity.  For  hygroscopic 
aerosol,  these  leqxmses  are  different  for  increasing  and  decreasing  relative  humidity, 
resulting  in  a  hysteresis  in  the  growth  curves  (e.g.,  by  Winkler,  1971).  On  the  other  hand,  a 
similar  behaviour  has  not  been  observed  in  die  relation  between  the  extinction  coefficient 
and  the  relative  humidity. 

In  many  cases,  die  visibility  conelates  bodi  with  die  badcscatter  and  with  the  extinction 
coefficients.  However,  dieie  are  also  periods  where  the  visibility  does  not  correlate  with 
either  of  diese  parameters.  This  depends  peihiqis  on  the  type  of  aerosoL 

Tbe  reference  values  for  the  badcscatter-  and  die  extinction  coefficients  are  calculated  by 
linear  regiessicHi  cm  horizontal  measurements.  Calculated  standard  deviations  in  tbe 
extinction  show  a  minimum  at  about  0.8  km-’,  and  for  die  badcscatter  at  about  0.02  km-i, 
vdien  die  range  of  die  lidar  is  at  niMcimiini.  A  mote  extensive  study  on  this  subject,  in  which 
die  quality  of  die  slope  medicxl  has  been  anatysed  as  a  funcdcm  of  the  signal-to-noise  ratio 
aod  die  actoal  extinction  will  be  publisbed  elaewbeie  (Kunz  and  DeLeeuw,  1992). 

Analysis  of  a  selected  groiqi  (no  precqiitatioo)  of  profiles  show  ditt  below  4(X)  m  die 
averaged  extinction  has  die  tendency  to  peak  d  about  02  km'>.  At  higher  altitudes  the 
variation  in  die  extincticn  increases,  while  cn  average  die  extinction  coefficient  decreases 

with  height 

Recommendations  fw  fnture  wodc: 

storage  of  die  raw  lidar  data  is  lecommoided  (aod  is  possible  nowadays  e.g.  with  large 
cspadty  storage  media)  to  nqirocess  die  data  widi  ahemttive  inversion  algorithms. 
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die  inb<nnogeoeity  of  the  bonzootal  extincdoa-  and  backscatter  coeffideots  can  be  expressed 
in  the  standard  deviations  as  calculated  widi  die  skpe  mediod. 

it  is  desirable  to  extend  the  altitude  range  to  2000  m  or  more  to  increase  the  probability  of 
cloud  detection. 

a  system,  bke  a  rain  gauge,  is  required  to  indicate  periods  with  prec^itation 

a  sinqile  photodetector  or  a  commercial  pyranometer  can  be  applied  to  measure  die  amount 
of  background  radiation  (and  periods  of  sun  shine)  for  estimating  die  amount  of  background 
noise  and  for  indicating  night  time  periods 

the  quality  of  logarithmic  ariqilifiers  should  be  investigated  in  more  detail  as  regards 
bandwidth,  noise  and  the  transfer  for  small  signals 

the  dynamic  range  and  the  noise  {noperties  of  avalanche  photodiodes  under  different 
background  conditions,  including  the  iq^lied  anylifiers,  should  be  investigated  to  estimate 
the  perfotmance  of  the  optical  receiver  (Kunz  and  Moerman,  1992) 

the  recording  frequency  of  the  meteo  station  (4  per  hour  during  die  project)  should  be 
adiqited  to  the  frequency  of  the  lidar  (10  per  hour) 

the  meteo  data  should  not  only  provide  the  mean  value  per  parameter  over  the  time  interval, 
but  also  the  standard  deviatirm  to  determine  the  variation  (implemented  at  this  moment) 
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APPENDIX  A;  RANGE  OF  UDAR  SYSmiS 

In  this  report,  the  range  of  a  lidar  syston  is  defined  as  that  range  where  the  received  power  is 
equivalent  to  the  noise  equivakot  power  of  the  receiver.  Noise  is  generated  by  the  eJectroaics  and 
by  the  background  radiation.  In  the  general  case  of  an  inhomogeneous  atmoqdiere,  it  is  not 
possible  to  predict  the  maximum  detection  range  but  in  die  qiecial  case  of  a  homogeneous 
atmosfdiere,  this  maximum  range  can  be  calculated  by  assuming  a  linear  relation  between  die 
backscatter-  and  extinction  coefficients.  The  starting  point  for  the  derivation  of  die  maximiim 
range  is  the  common  model  for  single  scatter  lidar  signals  fiom  a  hranogeneous  atmosidiere. 
Substituti<Mi  of  o((R)s  ot,  p(R>s  ^  and  the  noise  equivalent  power,  P^,  for  die  received  power  in 
equatkm  (3. 1)  leads  to: 


P.=E,-4 


2  4-rr  R^m 


‘opt 


^-2aR 

■e  m 


(A.l) 


where: 

A 

c 

E» 

P. 

R» 

Top. 

a 

P 


area  of  the  receiver  in  km^ 
veloc.  of  lightinkm/s 
laser  energy  in  J 

effective  noise  power  of  die  detector  in  W 
maximum  range  in  km 
transmission  of  die  optics 
spatial  extinction  coefficient  in  km'^ 
backscatter  coefficient  in  km** 


Taking  the  logaridun  of  (A.1)  leads  to: 


2.aR„+2.1„(R„)  =  l.(K,'.p) 


(A -2) 


where  is  the  lidar  system  constant  according  to: 
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I 

‘  "  2  4  il  P. 


(A -3) 


For  a  Rayleigh  atmosidtete  at  a  wavelengdi  1060  nm,  a  =  0.87- 10*^  km**  and  P  =  1.30*  10*^ 
km'^  SubstitutKMi  of  these  values  in  (AJ)  yields  a  maximum  range  for  the  Mini-Iidar  of  4S0  m. 

The  maximum  range  in  a  homogeneous  atmosphere  is  a  fimction  of  the  extinction  and 
the  backscatter.  If  a  ratio  C,  =  (generally  with  ksl)  is  assumed,  R„,..  can  be  ^aimiatBri  by 
taking  the  total  derivative  of  (A.2).  This  results  in: 


2*adR^+2R^da+*-^dR^—  (A*4) 

I'm**  ® 

The  extinction,  at  which  the  maximum  range  occurs,  is  found  by  explicitly  writing  dR^/da  and 
equating  this  term  to  zero.  This  results  in: 


(A  *5) 


Knally,  the  maximum  range  R^  is  found  by  substituting  (A.5)  in  (A.2) 
where  e  is  the  base  of  die  natural  logarithm. 


Equation  (A.6)  diows  that  an  increase  of  the  qrstem  power  by  flnee  orders  of  m«gtiitiirf<»  is 
required  to  inqnove  die  range  by  onty  one  order  of  magnitude. 
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Range  Uic  Mini-Uriif ; 

The  propoties  of  die  Mmi-Lidar  (after  a  veiy  tong  tune  of  operatkn)  are: 
A  »  4.310^kni2 

c  s  3.0  lO’  km/s 

a  SOmJ 
P,  =  8.3 10^  W 

Top.  -  0.5 

Substitation  of  dieae  values  in  (A.3)  provides  a  value  for  K': 

K'*  154.6  ktn^ 


With  C.S  0.07  and  facmula  (A.4)  the  maximum  range  is  found  to  be  1.26  km.  This  range  occurs 
when  die  horizontal  extinctioc  is  0.40 1an'‘  as  diown  by  formula  (A.S). 
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APFENDKB:  CALIBRATION  OF  THE  LOG  UNIT 


The  transfer  of  die  log  tmit  can  be  deacribed 


Uo-L,  +  L2ln(U,/U,);  (B.l) 

for  sufficient  large  ir^Nit  signals  U|, 
where: 

L]  «  constant  in  V 

L2  s  Constantin  V 

Up  -  ouqnit  voltage  in  V 
Uj  s  iiqmt  voltage  in  V 

U,  a  reference  voltage,  1  V 

The  reference  voltage  U,  is  added  to  make  in  the  argument  of  the  log  functkn  dimensionless. 

The  purpose  of  the  calibration  is  to  determine  the  constants  Lj  and  L2  and  their  accuracies  in  die 
interval  if^iere  die  transfer  is  logarithmic.  Therefcne,  a  square  wave  test  pulse  with  an  amplitude 
of  4.78  V  has  been  qiplied  at  die  input  of  die  amplifier.  A  passive  attenuator  with  a  range  from  0- 
1(X)  dB  (20  dB  is  a  factor  10  in  voltage),  in  stqis  of  3  dB,  is  used  to  cover  the  v^ok  irqnit  range. 
A  seamd  test  pulse,  with  an  anqilitode  of  20  %  of  die  main  test  pulse  and  smaller  in  width,  was 
added  to  investigate  the  dynamic  bdutviour  of  the  transfer  function.  MoteoveT,  diis  provided  a 
second  mediod  to  calculate  die  transfer  function. 
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Figure  B.l:  Reqxxite  of  the  log-unit  es  a  fiaictioo  of  the  eOenuator  (cttings.  The  sh^  of  the  test 
pulse  is  indicated  in  die  figure.  The  lower  set  of  data  punts  are  fiom  the  test  pul»  with 
an^rihude  1';  die  upper  data  set  nte  from  aio{ditude  7'. 

In  die  range  fiom  0  to  80  dB,  the  transfer  is  assumed  to  be  logarithmic;  in  dte  range  from  80  to 
100  dB,  the  transfer  is  s^oximated  by  sectmd  order  curve  fitting.  (Later,  it  appeared  that  for 
small  signals  the  transfer  becomes  linear.) 

The  constants  Li  and  detennined  by  linear  regressioo  analysis  over  the  range  0  to  80  dB 
(25  measuiemesit  points).  This  provides: 


L,  =377.29±  1.40V  and  (B.2a) 

Lj  =  44.73  ±0.20V 


(B.2b) 
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The  value  of  the  rawafnt  Lj  also  be  invoted  from  die  additifwuil  test  pulse.  This  cao  be  seen 
by  e:q>aiidmg  (B.l)  to: 


With  (B.l)  diis  results  in: 


L2  = 


AUs 


The  results  are  plotted  in  Figure  B.2. 


(B.3) 


(B.4) 
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FiguieB.2:  Calcu]atedvalueofL2*‘’<’<’'’‘^^^<l°**>**(^-^)' 

Analysis  of  die  data  over  die  input  range  from  0  to  80  dB  yields: 


(B.5) 


L,- 42.06  ±  4.19  V 
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It  is  clear  diat  diis  mediod  provides  a  similar  value  for  L2  as  in  (B^b)  but  with  a  larger  standard 
deviatkm. 

accuracvoftiieiiivertedinputsignals 

The  uncertainty  in  the  values  of  and  Lj  results  in  an  inaccuracy  of  foe  inverted  n^ut  signal 
fitom  Uo-  If  it  is  assumed  foat  foe  output  voltage,  U^,  can  be  measured  very  accnratety,  than  foe 
uncertainty  in  the  input  voltage  can  be  written  as: 


(B6) 

cm2 

EvaluaticRi  of  foe  partial  differentials  with  (B.l)  results  in: 


U,  (U.-L,)u, 

dUi  — - — dLi - - <11-2 

i-2 


which  leads  to: 


(B-7) 


dUj  dL>i  <11-2 

"uj  T7  T7 

Substitutioa  of  foe  numerical  values  of  foe  variables  leads  to  the  following  maximum  enors; 


ln(U,) 


(B-8) 


1.40  0.20 

44.75  44.75 


•ln(0.0003)  *0.067 


for  an  ii^t  voltage  of  0.0003  V:  and 
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dUi 

Uj 


1.40  ^  0.20 
44.75  44.75 


0.036 


-  for  an  iiqmt  voltage  of  3  V 


tran^  for  very  low  input  Mgnala 

Hgure  B.l  shows  Out  for  very  snull  input  aignal*  (snuller  than  5-10^  V)  die  transfer  is  no  longer 
logaridunic.  Therefore  a  second  order  curve  has  been  fitted  to  the  points  from  80  dB  to  100  dB. 
This  results  in  the  following  fit 


Uj  = -0.705397  +  5.34466.  U„  -6.242335.  Uo^:  (B.9) 

NOTE:  The  analysis  thus  far  is  based  on  a  transfer  of  the  log-unit  for  positive  signals  only, 
because  negative  input  signals  would  provide  unreal  outputs.  Later  measurements  showed 
however,  that  the  transfer  of  a  logaiithmic  amplifier  is  finite  for  very  small  signals.  This 
means  that  a  different  approach  is  necessary  as  explained  in  Kunz,  1990- 


Electrical  bandwidfli  of  the  log-unit 

The  transfer  function  of  the  logarithmic  aoqtlifier  as  a  function  of  frequency  can  in  fact  only  be 
justified  for  very  small  iiqnit  signals  ^^lere  tihe  transfer  is  linear.  Strtmger  hqmt  signals  produce 
hannonics  due  to  die  logarithmic  characteristics.  Neverdieless,  the  electrical  bandwiddi  has  been 
detennined  for  a  number  of  different  hqnt  kvds  and  has  been  compared  with  a  more  expensive 
log-unit  The  results,  which  are  shown  in  Figure  B3,  indicatB  that  die  bandwiddi  (3  dB  point)  is 
about  6  MHz  for  small  signals  and  increases  slight^  to  8  MHz  for  large  iiqiot  signals.  The  slope 
is  about  12  dB  per  octave  which  means  diat  signals  widi  high  frequency  cooqxments  are  strongly 
attenuated.  (A  log-unit  from  a  different  supplier  indicated  a  bandwidth  of  better  than  20  MHz 
widi  a  slope  of  less  than  2  dB  per  octave). 
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APPENDIX  C:  EFFECTS  OF  THE  TRIGGER  INACCURACY 

The  pupose  of  lidar  is,  among  others,  to  measure  die  extinction  and  the  backscatter  {rofiles  in 
the  atmosphere  as  a  function  of  range.  However,  die  signals  decrease  with  die  square  of  the  range. 
This  is  called  'geometric  attenuation'.  It  is  common  practice,  to  start  signal  processing  with 
compensating  for  this  geometric  attenuation.  The  moment  that  the  laser  pulse  is  emitted  is  the 
reference  for  die  starting  point  of  this  process.  However,  if  no  qiedal  precautions  are  takm,  there 
is  always  an  uncertainty  of  maximum  fins  sample  intovaL  Ftom  the  viewpoint  of  signal 
processing,  the  signal  is  corrected  with  (R±AR)2  in  stead  of  with  R^.  The  influence  of  this  error 
on  die  calculated  extinction-  and  backscatter  coefficients  is  discussed  in  diis  secticm. 

Because  it  is  not  possible  to  treat  all  possible  inhomogeneous  situations,  imly  signals  from 
homogeneous  atmosidieres  are  cmisideied.  Furthermore,  this  analysis  should  only  be  considered 
as  a  first  step  to  quantify  this  problem. 

Fot  simplicity,  it  is  assumed  that  the  lidar  signal  can  be  described  with  the  following  simplified 
model: 


Pd(R)=— 


-2«R 


where: 


K|  =  lidar  system  ctmstant  in  W-km^ 
Pj(R)  =  received  power  in  W 

R  =  range  inkm 

a  =  extinction  coefficfont  in  km'* 

P  =  backscatter  coefficient  in  km'* 


(Cl) 


If  there  is  a  shift  in  the  trigger  moment,  than  the  range  cmrection  is  performed  with  an  additional 
unknown  shift  AR  Ihis  can  be  described  as: 
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S(R)  =  P(R)•(R±AR)^  (C-2) 

where  S(R)  is  the  range  corrected  lidar  signal  in  W-km^. 

If  AR=0,  than  both  the  extincticMi  and  the  backscatter  can  be  determined  firom  a  linear  regression 
fit  on  LN[S(R)].  However,  if  AR  is  not  equal  to  zero,  the  {qtplication  of  linear  regression  to 
determine  the  errors  in  the  calculated  extinctioa-  and  backscatter  coefficients  leads  to  very 
complicated  analysis.  Therefore,  the  problem  is  split  into  two  parts: 

-  determine  the  error  in  the  backscatter  assuming  that  the  extinction  is  known 

-  determine  the  error  in  the  extinction  assuming  that  the  backscatter  is  known 


the  inverted  badasater 

The  backscatter  coefficient  is  determined  by  averaging  over  interval  Rj  to  Rj  according  to: 


P±AP  = 


Rj 

J 


Pd(R)(R±AR)^ 


(Rj-Rj)  J  Ki-e 


-2.«R 


dR 


(C-3) 


For  small  values  of  AR/R,  the  quadratic  term  is  approximated  by  the  first  term  of  the  series 
expansion.  Than  the  integral  leads  to: 


Ap  2  AR 

P  (R2  —  Ri) 


(C-4) 


Thus  the  errOT  in  the  backscatter  is  proportional  to  die  trigger  inaccuracy  and  dq>ends  also  on  the 
position  and  the  length  of  the  selected  interval. 
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For  exanq)le,  a  tciggex  inaccuracy  of  30  ns  leads  to  an  enor  of  about  1.2  %  in  tte  badcscatter 
(RjslOOm  and  Rj^lOOOm).  The  error  increases  to  about  3.1  %ifR2is200m. 

Error  in  the  inverted  extinction 

The  error  in  die  extinction  is  determined  by  conjuring  tbe  average  slope  of  the  selected  interval 
with  the  actual  extinction  coefficioit  Furthermore,  it  is  assumed  that  the  backscatter  P  is  knowa 
Than  one  can  write: 


a±  Aa= 


1 

(R2  ~Ri) 


7  -  u. 

■pjfRlAR)^' 

J2  dR 

PK, 

•dR 


(C-5) 


For  small  values  of  AR/R,  the  quadratic  term  can  be  tqiproximated  by  the  first  term  of  the  series 
expansion.  Subsequent  evaluation  of  tbe  integral  leads  to  tbe  following  relative  error  in  the 
extinction: 


^ - ^  (C.6) 

a  aRiR2 

For  example,  if  the  trigger  inaccuracy  is  30  ns  and  the  measuranent  is  made  in  a  clear 
atmosidiere  with  0;=  0.1  km'^  the  error  in  the  extinction  becomes  56  %  when  Ri=  100  m  and  R2= 
800  m! 


TNOraport 


AppeodixD 


Page 

D.l 


APPENDIX  D:  EFFECTS  OF  THE  TRANSFER  ACCURACY  OF  THE  LOG  UNIT 

The  uncertainty  of  the  transfisr  of  the  log  unit  leads  to  etiors  in  the  calculated  extinction-  and 
backscatter  coefficients.  For  sinqilicity,  (mly  signals  from  hmnogeneous  atmosiriiaes  are 
considered.  The  problem  is  ^lit  in  two  parts: 

/•ainiiatf.  the  influence  on  die  extinction,  assuming  that  die  backscatter  is  known 
calculate  the  influence  on  the  backscatter,  assuming  that  the  extinction  is  known 

The  analysis  starts  with  the  sin^lified  model  of  the  lidar  equation  from  homogeneous 
atmospheres: 


U,(R)  = 


K.D-P 

R^ 


.g-2-aR 


(Dl) 


were: 

D 

sensitivity  of  the  detector  in  V/W 

K,  = 

lidar  system  ctmstant  in  Wkm^ 

R 

range  in  km 

a  = 

extinction  coefficient  in  km** 

P 

backscatter  coefficient  in  km** 

According  to  equation  (B.l),  the  output  of  the  log  unit  becomes: 


Uo(R)  =  Li  +L2  [LNIK,  .D  p)-2  o  R-2.LN(R)] 


(D.2) 


Note  diat  the  referaice  voltage  has  been  set  to  U,®  1  V  and  will  therefore  not  qipear  in  this 
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equadoa 

The  Mtmctifm  coefiEicient  a  can  explicitly  be  written  as: 

LN(K)i  •  D  •  p)  -  2  •  LN(R)  -  ^ 

a  = - — - h -  (D.3) 

2R 

Now  the  total  eiror  in  die  extinction,  due  to  errors  in  die  transfer  (xxistants  of  the  log  unit,  can  be 
written  as: 


,  3a  3a  „ 


with: 


3a _ 1 

3Li  *2  R-L2 


3a  _  LN(Ui) 
3L2  2RL2 


(D.4) 


(D.5a) 


(D5h) 


The  total  error  in  the  exdnctimi  becomes: 


2  R  Lj  ‘  2  R  L 


(D.6) 


The  average  error  in  the  extinction  over  interval  Rj  to  R2  becomes: 


Aa= 


(R: 


1  y  AL,  ^  LN(U0-AL,1 

^  iL2  R  L2  2  R.L2  J 


dR 


(D.7) 


After  evaluating  eiqnession  (D.7)  one  Snds; 


Aa  = 


AL,LN(^)  +  AL2 
Ri 


LN(Ki  •  D .  p).  LN(^)  -  2.|^LN(^) J  -  2  •  a  •  (Rj  -  Rj ) 


2  Lj  (Rj-Ri) 


(D.8) 


The  relative  uncertainty  in  o  h«.<  been  calculated  by  substituting  the  values  of  K],  D,  Li,  AL],  1^, 
ALj  and  file  relation  of  p/o  in  (D.8).  (See  Table  2.1  and  Appendix  B;  for  p/a  a  value  of  O.Cn  has 
been  substituted).  For  R,  a  value  of  0.1  km  has  been  chosen  and  for  Rj  a  value  of  1  km.  The 
relative  uncertainty  in  the  extinction,  Aajet,  due  to  uncertainty  in  L|  and  1^  amonnts  -385 
%,  -23.1  %  and  -1.33  %,  for  a  values  of  0.01  km‘^  0.1  km'*  and  1.0  km**  respectively.  These 
ca^n1lllrion«  show  that  it  is  difficult  to  measure  small  extinction  values  with  diis  log  unit  (and  this 
stnaii  lidar  system).  Moderate  and  large  extinctions  can  be  measured  more  accurately.  Note  that 
the  uncertainty  can  be  even  worse  if  the  detector  noise,  the  spatial  variation  of  the  annosphcric 
scattering  and  the  uncertainty  in  the  other  parameters  are  taken  into  account 

In  the  way  the  error  in  the  backscatter  aieffident  is  calculated,  assuming  fliat  the  extinction 
coefficient  is  known. 

with: 

^  =  (D.lOa) 

aLi  Lj 

and 

8p  -LN(Ui)  (DlOb) 

9Li2  I.<2 

The  total  uncertainty  in  the  backscatter  coefficioit  becomes: 

Lj  L2 


The  averaged  enm  in  die  backscatter  over  interval  Rj  to  R2  becomes: 
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AP  = 


-P 

(Ra  -Ri) 


.  LN(U,) 

+  — - - ALj 


(D.12) 


After  evaluating  expressitm  (D.12)  the  absolute  uncertainty  in  the  backscatter  beccunes: 


(D.13) 


AP  =  -P 


AL, 


+  ^jLN(KiDp)-a(R2+Ri)-2-2-^-2-^-^^?^  RiLN(R,) 


(R2-R1) 


The  relative  uncertainty  in  the  backscatter  has  been  calculated  with  the  same  system  values  as 
described  under  (D.8).  This  results  in  an  uncertainty  of  1.32  %  in  Ap/p  if  the  extiiKtion  is  0.01 
km'*.  For  extinctions  of  respectively  0.1  and  1.0  km'*  the  relative  uncertainty  is  1.9  %  and  2.9  %. 
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Qnclusion 

The  uncertainty  in  the  tiansfin  function  of  the  iqjplkd  log-anq>]ifier  has  only  a  sli^  influence  on 
the  calculated  backscatter.  The  influence  on  Ae  calculated  extinction  is  much  stronger  and 
increases  as  the  extinction  decreases.  Unoettainties  as  large  as  a  factor  3  are  possible  if  the 
extinctkm  is  0.01  km'^ 
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APPENDIX  E:  DEFINITION  OF  THE  PARAMETERS  IN  THE  DATA  BASE 


The  following  paiameten  are  pfcseat  in  die  final  dau  base: 

1  -  lidardate 

2  -  meteodniB 

3  -  lidartime 

4  -  number  of  lidar  reference  measnrements 

5  -  meteodate 

6  -  sampling  interval  lidar 

7  -  meteocode 

8  - 

9  -  cloud  detection  on/off 

10  -  tenyierature 

11  -  dewpoint 

12  -  relative  humidity 

13  •  lidar  elevation  angle 

14  -  pressure 

15  •  wind  velocity 

16  •  wind  direction 

17  •  visibility  AEG 

18  -  radon  counter 

19  -  transmission  over  265  m 

20  -  minimum  value  profile 

21  -  maximum  value  profile 

22  -  min,  altitude  profile 

23  •  max.  altitude  profile 

24  laser  energy 

25  -  horizontal  extinction  in  km* 

26  horizontal  backscatter  in  km'> 

27  -  horizontal  delta  extinction  in  km~* 

28  -  horizontal  delta  backscatter  in  knr* 

29  -  turbulence 

30  -  cloud  cover 

31  -  running  averaged  exdnctioa 

32  -  running  averaged  backscatter 

33  -  altitude  reliability?? 

34  -  cloud  altitude  in  m 

35  -  cloud  thickness  in  m 
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APPENDIX  F:  SOME  RESULTS  OF  THE  SODAR  SYSTEM 

In  1984  and  1985,  an  experimental  sodar  system  has  been  constnicted  widi  a  high-power  q)eaker 
and  a  parabolic  disk  widi  a  diameter  of  about  60  cm.  The  transmitter  openued  at  a  frequency  of 
1650  Hz  widi  an  acoustic  power  of  50  W.  The  sodar  syston  was  mounted  on  the  roof  of  the 
institute  at  an  altitude  of  about  15  m  above  ground  leveL  A  small  wooden  cylinder  around  die 
system  should  stqipress  the  background  noise.  The  bandwiddi  of  the  receiver  was  about  40  Hz 
which  was  realized  widi  a  boxcar  integrator.  A  PET  2000  Commodore  Business  Conqiuter 
controlled  the  system  and  the  data  were  stored  on  audio  cassette  tape  which  had  a  storage  capacity 
of  about  100  signals.  Therefore  each  sodar  profile  presented  here  is  an  average  of  about  15 
minutes. 

Some  of  the  obtained  sodar  data  are  presented  in  Figure  F.l.  Here  the  received  power  has  been 
coded  in  false  color  in  a  time  versus  height  diagram.  The  vertical  axis  represents  the  height  from 
15  to  415  m  whereas  die  time  has  been  plotted  along  the  horizontal  axis.  The  hourly  intervals 
have  been  indicated  by  the  small  vertical  tics  at  the  bottinn  and  the  tiqi  of  the  figure. 

Two  different  tune  intervals  can  be  distinguished  viz.  from  07:00  am.  to  07:00  pm.  and  from 
07:00  pm  to  07:00  am  As  in  many  odier  data  files,  the  data  show  much  structure  during  the  day 
from  about  07:00  a.m.,  when  die  laboratory  opens,  until  06:00  pm.  when  the  labaratoty  closes.  It 
is  conjectured  dierefore,  diat  the  measured  structure  is  not  only  caused  by  diermal  effects  but  also 
by  the  background  noise  of,  e.g.,  die  engines  of  cars  and  other  human  activities.  During  the 
evening  and  the  night,  die  background  noise  is  at  a  lower  level  and  some  horizontal  structure 
becomes  visible  between  about  50  m  and  300  m  altitude. 

The  results  of  a  two  diys  period,  which  are  shown  in  Hgure  F.lb,  show  similar  effects.  The 
horizoiital  layering  is  clearty  visible  during  die  night  but  disqipears  here  after  07:00  am 

The  sodar  results  do  not  show  any  coirelatkm  widi  die  lidar  results.  This  can  partly  be  eiqilaiiied 
fi'om  the  different  interaction  processes  and  die  limited  range  of  the  sodar  system  For  possible 
future  sodar  work,  it  is  recommended  to  have  abetter  acoustic  isolation  from  the  mvironment  and 
a  imne  sensitive  system 
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Figure  F.l:  Veitical  (tnicture  of  the  auaospben,  coded  in  false  odor,  as  measured  widt  an  acoustic 
sounder  (sodar).  Veitical  axis  0-400  m;  horizontal  axis  is  time:  the  smaD  vertical  lines 
indicate  Ore  houily  mterval  and  the  laige  veitical  lines  indicate  midni^t  Sodar  results  from 
FUvuaiy  4  to  8 ,1^.  (a)  and  sodar  results  from  Feteuaiy  11  to  12, 1985  (b). 
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